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£vq ABSTRACT 

Context. Large spectroscopic surveys have discovered very peculiar and hitherto unknown types of active galactic nuclei (AGN). 

^vj Such rare objects may hold clues to the accretion history of the supermassive black holes at the centres of galaxies. 

Aims. We aim to create a sizeable sample of unusual quasars from the unprecedented spectroscopic database of the Sloan Digital Sky 

Survey (SDSS). 

Methods. We exploit the spectral archive of the SDSS Data Release 7 to select unusual quasar spectra. The selection method is based 

^h on a combination of the power of self-organising maps and the visual inspection of a huge number of spectra. Self-organising maps 

^. were applied to nearly 10 5 spectra classified as quasars at redshifts from z = 0.6 to 4.3 by the SDSS pipeline. Particular attention was 

paid to minimise possible contamination by rare peculiar stellar spectral types. All selected quasar spectra were individually studied 
to determine the object type and the redshift. 

i—, Results. We present a catalogue of 1005 quasars with unusual spectra. These spectra are dominated by either broad absorption lines 

(BALs; 42%), unusual red continua (27%), weak emission lines (18%), or conspicuously strong optical and/or UV iron emission 

_|_J (11%). This large sample provides a useful resource for both studying properties and relations of/between different types of unusual 

quasars and selecting particularly interesting objects, even though the compilation is not aimed at completeness in a quantifiable sense. 
The spectra are grouped into six types for which composite spectra are constructed and mean properties are computed. Remarkably, 
all these types turn out to be on average more luminous than comparison samples of normal quasars after a statistical correction is 
made for intrinsic reddening (E(B-V) « to 0.4 for SMC-like extinction). Both the unusual BAL quasars and the strong iron emitters 
have significantly lower radio luminosities than normal quasars. We also confirm that strong BALs avoid the most radio-luminous 
quasars. For 32 particularly interesting objects, individual spectra are presented. Among these objects are quasars with many narrow 
BAL troughs and one quasar where the continuum is strongly suppressed by overlapping BAL troughs across nearly the whole SDSS 
spectrum. Finally, we create a sample of quasars similar to the two "mysterious" objects discovered by Hall et al. (2002) and briefly 

l — — ' discuss the quasar properties and possible explanations of their highly peculiar spectra. 
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1 1. Introduction UV-optical SEP is well represented by global co mpos ite spectra 

S— : produced from large samples (e.g., Francis et al. 1991 ; Zheng et 

v-2 From the very beginning ofquasar astronomy (Schmidt [1963J al [J9921 Vanden Berk et al. 2001 >, where discretion should be 

CO Greenstein & Matthews [1963| Sandage LL965J, optical spectra exercised however when using the composite of one sample as 

O have been one of the most important sources of information a tem pi ate for other samples (Vanden Berk et al. 2001 ; Richards 

CN about this class of fascinating objects. Quasars are the most lu- et a j 2003 ) 
,— ' minous types of active galactic nuclei (AGNs). The average ul- 

> traviolet (UV) and optical spectral energy distribution (SED) of II is t0 be expected that various predicted spectral features, 

H AGNs is offen found to be widely consistent with the predic- such as the s P ectral slo P e or the stren g th of broad emission lines, 

X tions of the standard model of a geometrically thin and optically can chan § e from quasar to quasar owing to variations in physical 

H thick accretion disc around a massive or supermassive black hole Parameters (e.g., Koratkar & Blaes [1999J Laor & Davis ESD, 

(Shakura & Sunyaev 1973 ; Shields ElS Laor 1990; Blaes et and there is empirical evidence that spectral properties can also 

al. 2001 \ The continuum radiation field from the accretion disc be modified by variability (Kinney et al. 1990; Wilhite et al. 

provides the photons for the ionisation of the surrounding gas 2005 ; Hal1 et aLEHU Meusinger et al.ESD. Further modifica- 

and powers broad and narrow emission lines. The unified AGN tions of the spectrum can be caused by absorption from material 

model (Antonucci 1993 ; Urry & Padovani 1995 ) allows us to un- relate d to the quasar itself, to the host galaxy, or from unrelated 

derstand some differences in the spectra of various AGN types intervening matter along the line of sight. In addition, some spec- 

as being caused by different orientations of the sightline towards tral features can be_diluted by unrelated emission components 

the central engine. For the vast majority of known quasars, the (Johnston et al. [2003J Kishimoto et al. 2008) or magnified by 

gravitational lensing (Abajas et al. 2007; Yonehara et al. 2008; 

~~Z, ~ ; ZTT, ,,, . , -,,,-, ■ x- Blackburneetal. 201J). 

The lull catalogue (Table 3) is only available in electronic form at 

the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via On the other hand, it is not to be expected that quasars are 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-i-A/ fully described by a static model, because "quasars are not things 
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so much as processes'" (Richards et al. 120111) . Over approxi- 
mately the past decade, spectroscopic surveys have revealed a 
variety of quasar spectra that dramatically differ from the stan- 
dard SED, and confirmed the existence of populations of exotic 
and hitherto unknown quasar types (Fan et al. 1992 White et al. 
[20001 Menou et al I20UT1 Hines Egg D Hal l et al. [20021 [2004; 
Brunner et al. [20031 Meusinger et al. [20051 Plotkin et al. 2008; 
Diamond-Stanic et al. [2009; Urrutia et al. 120091) . The discov- 
ery and investigation of these rare peculiar objects is important 
because they may represent links to special evolutionary stages 
and provide an excellent laboratory to study feedback mecha- 
nisms between star formation and accretion activity. There are 
strong indications from different lines of thought for a tight rela- 
tionship between the activity and growth of massive black holes 
and the evolution of their host galaxies (e.g., Sanders et al. 1988 ; 
Magorrian et al. [T9981 Kau ffmann & Hahnelt 2000; Taniguchi 
I20U51 Ha ring & Rix[2004; Hopkins et al. 2006; Kauffmann & 
Heckman 2009 ; Meng et al. l20T0l Shabala et al. [20TT1) . Much at- 
tention has been focused in particular on identifying the youthful 
quasar population. 

There are several tentative arguments supporting the view 
that young quasars are associated with gaseous outflows and 
absorption by dust and gas. Many AGNs show blueshifted ab- 
sorption from outflowing matter (Ganguly & Brotherton 120081) . 
Very broad absorption lines (BALs) corresponding to a wide va- 
riety of velocities are observed in luminous quasars. These BALs 
are naturally explained as footprints of powerful, sub-relativistic 
(up to ~ 0.2c) outflows (Weymann et al. 1991), which are per- 
haps related to radiation-driven winds from the accretion disc 
(Murray & Chiang 1998). Broad absorption lines are thus the 
most obvious manifestation of matter accelerated by the cen- 
tral engine of AGNs. It has been frequently argued that the 
BAL fraction can be explained as being caused by an orienta- 
tion effect within the unification model where otherwise normal 
quasars are seen at sightlines close to the equatorial plane of 
the AGN (Weymann et al. 119911) . Typical BAL quasars are often 
subdivided into those showing broad absorption lines from only 
high-ionisation species such as C iv and N v (HiBALs) and those 
showing in addition also absorption from low-ionisation species 
such as Mg ii and Al m (LoBALs). A rare subclass of the LoBAL 
quasars are the FeLoBALs with absorption from metastable ex- 
cited states of Fell and Felll (Hazard et al. 119871 Becker et 
al. 12000b . The BAL quasars are thought to account for between 
10% and- 40% of all qua sars (H ewett & Foltz 2003; Reichard 
et al 120031 Carball o et al . 120061 Tru mp et a l. [20061 Dai et al. 
120081 Gibson et al. [20091 Allen et al. [20TTT) . We note that the 
exact BAL fraction is not trivial to determine because of dif- 
ferential selection effects between BAL quasars and non-BAL 
quasars. In addition, different definitions of BAL quasars have 
been used in different studies. According to the criteria applied 
by Trump et al. ( 2006 ), BAL quasars comprise 26% of quasars, 
LoBALs about 1.3%, and FeLoBALs about 0.3% (these being 
the raw fractions without correction for the additional reddening 
of BAL quasars compared to non-BALs). 

Outflows are believed to be fundamental to understanding 
both the overall picture of AGNs and the AGN feedback. Despite 
their importance, BAL quasars remain poorly understood. In par- 
ticular, there is currently no concensus about the role of evo- 
lutionary effects in determining the properties of BAL quasars. 
The properties of HiBAL quasars appear to be consistent with 
the pure orientation model (e.g., Surdej & Hutsemekers 1987; 
Gallagher et al. [20071 Doi et al. [20091) . On the other hand, 
LoBAL quasars may not be explained by orientation alone (e.g., 
Ghosh & Punsly 2007 ; Montenegro-Montes et al. 12009 ; Zhang et 



al. l20T0l but see Gallagher et al. [20071 and Hall & Chaiet l20TTl) . 
Their properties probably depend on both orientation and evolu- 
tion effects (Richards et al. 201 1 ; Allen et al. 201 1 ). Based on the 
shape of the continuum in the UV and X-rays and on the BAL 
fraction in the infrared, it has been suggested that FeLoBAL 
quasars in particular probe a much more obscured quasar popu- 
lation than non-BAL quasars and may represent a young evo- 
lutionary stage when a thick shroud of gas and dust is being 
expelled from the central region of the AGN (Voit et al. 1993; 
Becker et al.[T997l Canalizo & Stockton l200Tl Lipari et al. [20091 
Dai et al. 20 ll Farrah et al. [20T0b . 

It was suggested by early studies that AGNs spend a sub- 
stantial fraction of their lifetimes in a dust-enshrouded environ- 
ment. This picture was based on the scenario of merger-driven 
AGN activity (e.g., Sanders et al. 1988; Canalizo & Stockton 
ED QII Ko mossa et al.[2003; Hopkins et al. 12005120061 Guyon et 
al. 2006; Bennert et al. 2008 ). Major mergers of gas-rich galax- 
ies are very efficient in driving gas and dust to the central re- 
gions thus feeding both circumnuclear starbursts and accretion 
onto massive black holes (e.g., Hernquist 1989; Springel et al. 
2005 ). The natural consequences are not only a connection be- 
tween nuclear accretion activity and starbursts but also the pre- 
diction of a substantial fraction of very dusty quasars. Studies 
in the radio, infrared, and hard X-ray domain (e.g. Webster et 
al. H9 95; Cutri et al. 2001; White et al [20031 Glikman et al. 
2007; Leipski et al. [20081 Polletta et al. 2008) have revealed a 
large population (perhaps > 50%) of dust-reddened quasars that 
optical colour surveys tend to miss and that are not fully ex- 
plained by obscuration in the framework of the (static) unified 
model (Martinez- Sansigre et al. 2006; Georgakakis et al. 2009; 
Urrutia et al. 2009 ). Though it is much more difficult to find dust- 
obscured quasars in the optical, it is also not unrealistic to expect 
moderately reddened quasars in optical surveys. 

Another category of unusual quasar spectra that has been 
discovered by large spectroscopic surveys, typically display the 
unobscured UV quasar continuum but extremely weak or un- 
dectable emission lines. The first such quasars were detected 
more than one decade ago (McDowell et al. 1995; Fan et al. 
1999). About 80 luminous sources of this type with 2.2 < z < 5.9 
have since been identified by the Sloan Digital Sky Survey 
(Shemmer et al. 120101 and references therein). Various sugges- 
tions have been made to explain the absence of broad lines in 
unobscured quasars (e.g., Shemmer et al. I2010[ Laor & Davis 
2011 ; and references therein) but this AGN type remains puzz- 
ling. A subset of quasars with weak UV emission lines have been 
identified that are known as PHL 1811 analogs (Leighly et al. 
2001 ; Leighly et al. [20071 Wu et al. 1201112012b . Leighly et al. 
(2007 ) suggested that the weak UV emission lines of these ob- 
jects are due to an unusual soft SED deficient in ionising pho- 
tons, perhaps due t o ex ceptionally high Eddington ratios (see 
e.g., Shemmer et al. 2009 and references therein). On the basis 
of the multi-epoch photometric data of the quasars in the SDSS 
stripe 82, Meusinger et al. (120111) found that weak-line quasars 
are less variable than normal. This seems to be consistent with 
the suggestion (Hryniewicz et al. 2010) that these quasars are in 
an early evolutionary stage. However, we note that the sample 
of weak-line quasars with reliable variability data is small and 
needs to be checked carefully for possible stellar contamination. 

Unusually strong Fe n emission is observed for a small frac- 
tion of quasars. The Fen emission extends from the near-UV to 
the optical region and is often one of the most prominent spec- 
tral features. The strength of the optical Fen emission is one of 
the parameters that control a set of correlations between vari- 
ous emission line properties, known as the Eigenvector 1 (El, 
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Boroson & Green 1992). Eigenvector 1 is thought to ser ve as 
a surrogate H-R diagram for AGNs (Sulentic et al. 120071 ). The 
principal driver of El seems to be the accretion rate, which may 
be related to the evolutionary state. Strong iron emission is cor- 
related to the occurrence of BAL troughs (Boroson & Meyers 
[T951 Zhang etal. IMS. 

We finally note that double-peaked broad emission lines are 
another type of rare feature in quasar spectra (Halpern et al. 
I200T1 Strateva etal.l2003l Luo et al.|2009 ; Chornock et al. l20T0b . 
Strateva et al. (2003) presented a sample of 116 double-peaked 
B aimer line AGNs (corresponding to ~ 4% of low-z AGNs). The 
spectral signature of this subclass is displaced red and blue peaks 
of the emission lines similar to double-peaked emission lines 
found in the spectra of cataclysmic variables. Double-peaked 
emission lines are believed to originate from rotational motion 
in a relativistic accretion disc. High rotational velocities in com- 
bination with a highly inclined disc and strong reddening can 
alter the AGN spectrum significantly (Hall et al. 2002). To date, 
the double-peak effect has only been observed in low-ionisation 
lines, while the UV lines, which dominate the spectra of high- 
z quasars, are usually single-peaked. For some double-peaked 
AGNs, a subparsec supermassive black hole binary has been 
suggested as explanation (Gaskell 1983; Tang et al. 2009]). The 
existence of close massive black-hole binaries appears as a nat- 
ural consequence of galaxy mergers. 

About one decade ago, only a few peculiar quasars had been 
discovered, mainly from follow-up of radio sources from the 
FIRST survey (Becker et al. 11995120001 W hite et al. [20001) . The 
Sloan Digital Sky Survey (SDSS; York et al. l2Q06b demonstrated 
in its early phase that these rare objects are not unique, but 
are members of populations of very peculiar quasars that were 
occasionally found to have spectral properties never seen be- 
fore (Hall et al. 120021) . Meanwhile, the SDSS has dramatically 
increased the number of spectroscopically confirmed quasars. 
With its efficient quasar selection technique and huge number 
of high-quality spectra, the SDSS is uniquely qualified to signif- 
icantly increase the samples of unusual quasars. 

The present paper aims to select spectroscopic outliers 
among the quasars from the SDSS Seventh Data Release (DR7; 
Abazajian 2009). The database is shortly described in Sect.[2] 
The selection itself is based on an artificial neural network algo- 
rithm that uses unsupervised learning (Sect. [3]). The strategy of 
unusual quasar selection, the rejection of contaminating sources, 
and the classification of the selected objects are described in 
Sect.[4j The final catalogue and sample-averaged properties, in- 
cluding composite spectra, are discussed in Sect.|5j followed by 
the presentation and discussion of 32 particularly strange spec- 
tra. Summary and conclusions are given in Sect.[7] Standard cos- 
mological parameters Hq = 71 km s" 1 Mpc _1 ,Q m = 0.27, and 
Qa = 0.73 are used throughout the paper. 

2. Database: Sloan Digital Sky Survey 



The Sloan Digital Sky Survey (SDSS; York et al. [20061) mapped 
one quarter of the sky and performed a redshift survey of galax- 
ies, quasars, and stars. The total spectroscopic survey area covers 
9380 square degrees of the high Galactic latitude sky. The spec- 
tra were taken with the 2.5 m SDSS telescope at Apache Point 
Observatory equipped with a pair of double fibre-fed spectro- 
graphs. As the spectra are collected, the data is passed on and 
filtered through the spectroscopic pipeline. The automated sys- 
tem produces wavelength and flux-calibrated spectra, measures 
emission and absorption lines, classifies spectra, and determines 
the redshifts, which are, among other data, written into the head- 



ers of the fits files containing the spectra. The spectra are made 
publicly available via the SDSS Data Archive Server (DAS). 
The present study uses quasar spectra from the SDSS Seventh 
Data Release (DR7; Abazajan et al. 2009) available from the 
DAS. The DR7, marking the completion of the original goals 
of the SDSS and the end of the phase known as SDSS-II, con- 
tains 1 440961 spectra in total (after removing skies and dupli- 
cates), including 121 363 quasar spectra. The spectra cover the 
wavelength range 3 800-9200 A with a resolution of ~ 2000 and 
sampling of ~ 2.4 pixels per resolution element. 

The majority of quasars were selected for spectroscopic 
follow-up according to their deviation away from the "stellar lo- 
cus" in the multidimensional SDSS colour space. The SDSS fil- 
ter system allows us to discover quasars over the redshift range 
from ~ to 7. In addition to the multicolour selection, a small 
fraction of unresolved objects were selected as primary quasar 
candidates because their SDSS positions closely match either the 
positions of radio sources from the VLA FIRST survey (Becker 
et al. [T9951) or of ROSAT X-ray sources (Voges et al. 2000). 
A fully detailed description of the SDSS method of selecting 
quasars is given by Richards et al. (12002b . 

A catalogue of 105 783 quasars from the SDSS DR7 was 
compiled and published by Schneider et al. (120101) as the Fifth 
Edition of the SDSS Quasar Catalogue (hereafter QCDR7). This 
catalogue illustrates the unprecedented impact that SDSS has 
made for the quasar database. Nevertheless, since the present 
study is aimed at the search for very unusual SDSS quasar spec- 
tra, our quasar selection is based on the spectra available via 
the Data Archive Sever of the SDSS DR7 instead of using the 
QCDR7. 



3. Kohonen self-organising maps 



The self-organising map (SOM) algorithm (Kohonen 1982, 
2001 ) is an example of unsupervised learning of an artificial neu- 
ral network. It is an effective tool to project higher-dimensional 
input data onto a two-dimensional (2D) topological map. In this 
way, the algorithm has the capability of sorting input data onto 
a plane according to similarity and thus ideally prepares the data 
for visual inspection. To achieve this goal, each input element is 
assigned to one neuron in the network. All neurons of the net- 
work are organised onto a 2D grid. In each learning step, the in- 
put element is reassigned to some neuron, depending on match- 
ing criteria described in Sect. 3^2 Through this iterative stochas- 
tic process, the neural network is capable of learning, since it is 
designed to adapt to the input space. Therefore, input elements 
are not placed randomly, but increasingly closer to similar input 
elements, thereby clustering the input data. 

The characteristic evolution of the Kohonen map of quasar 
spectra is illustrated in Figs.[l]and[2] About 5 10 3 quasar spectra 
from the Fourth Data Release (DR4; Adelman-McCarthy et al. 
2006) were arranged in 200 iteration stepq^] The 18 panels show 
the 2D arrangement of the spectra at different iteration steps (top: 
1, 2, 3, 4, 6, 8, middle: 10, 20, 30, 40, 60, 80, bottom: 100, 120, 
140, 160, 180, 200). In Fig.[T] the redshifts z from the spectro- 
scopic pipeline of the SDSS are highlighted by means of colour 
coding. The gradient from dark to bright (black to yellow in the 
colour image) represents the range from low to high redshift. 
About 20% of the pixels are not associated with spectra and are 



1 The only reason for the use of DR4, instead of DR7, is that it can 
be analysed using substantially less computing time. As this map was 
produced for illustrative purposes only, we therefore chose to use the 
smaller database from the DR4. 
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Fig. 1. Temporal evolution (left to right then top to bottom) of the Kohonen SOM of 5 10 3 SDSS quasar spectra. The quasar redshifts 
are highlighted by means of colour coding. The gradient from dark (black) to bright (light yellow) represents the range from low to 
high redshift. 



marked in grey (see below). It is clearly seen that different red- 
shifts are separated in the final Kohonen map. 

A usual way of representing a SOM is the unified distance 
matrix (U-matrix; Ultsch & Siemon l 1990b shown in Fig.[2]for the 
same iteration steps as in Fig.[T] The U-matrix visualises the dif- 
ferences between the spectra and their neighbours. Light colours 
indicate high degrees of variation between adjacent spectra. 
Strong differences between the neighbours on the map are mea- 
sured, in particular, for rare peculiar spectra. Fig. [2] illustrates 
that these outliers tend to settle at the edges and corners of the 
final SOM. This is a useful property of the Kohonen maps that 
allows the efficient selection of unusual spectra. 

For the present study, the quasar spectra from the SDSS DR7 
were clustered according to their relative differences. At a given 
redshift, these differences are dominated by the shape of the con- 
tinuum and the presence of strong/broad absorption and emis- 
sion features. Ideally, the spectra will be sorted according to 
their spectral subtypes, thus differentiating between quasar sub- 
classes. 

In this Section we describe the preparation of the spec- 
tral dataset and the SOM algorithm in a mathematically hand- 
waving manner. A more detailed description of the method and 
the visualisation of results from the Kohonen mapping of large 
samples of spectra will be given elsewhere (in der Au et al., in 
preparation). 

3.1. Dataset 

An SQL query asking for quasars (spec_cln = 3 or 4, i.e., 
quasar or high-z quasar) with redshifts z > 0.5 in the SDSS 
DR7 returned 103 955 spectra. We are aware of the possibility 
that some unusual quasar spectra were classified by the SDSS 



spectroscopic pipeline as unknown objects (spectral class = un- 
known, spec_cln = 0). However, the subsample of the un- 
knowns is quite heterogeneous in nature and contains a high 
fraction of spectra with low signal-to-noise ratios (S/Ns). In ad- 
dition, there are, of course, no reliable redshifts given for these 
objects. As a consequence, analysing the SOMs of this sub- 
sample requires a different approach than for the quasars. The 
present paper is co ncerned with the SOMs of the quasar sample 
(but see Sect. |6.2| for an exception). The computation and sys- 
tematic analysis of the SOMs for the unknowns is in preparation 
and will be the subject of a separate investigation. 

As a first step, the overall size of the data set was reduced. 
We extracted the spectra from the fits files, together with the 
selected header keywords z , spec_cln, mjd, fiberid, and 
plate id. As our search for unusual quasars was designed to 
identify spectral peculiarities covering a wide spectral range, 
such as strong BALs, overlapping absorption troughs, and strong 
reddening, we reduced the resolution of the spectra used for the 
Kohonen method by a factor of four, which resulted in ~ 975 
pixels covering the wavelength range of 3800 - 9200 A in the ob- 
server frame. The consequence was a significant gain of compu- 
tation time. In the subsequent process of evaluating the selected 
spectra we used, however, the original SDSS spectra. All spectra 
were normalised to the integrated flux density. 

The Kohonen method applied to the whole sample of quasar 
spectra is expected to separate the spectra according to their red- 
shifts z (Fig. [I]). However, we repeat that it was our aim to search 
for unusual spectra. It is therefore advisable to apply the method 
to the quasars in narrow z intervals. The bin size Az should be 
small enough to ensure that the differences between the spec- 
tra, as seen by the SOM, caused by their different redshifts are 
smaller than the differences due to any spectral peculiarities. On 
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Fig. 2. As Fig.[lJ but showing the U-matrices. 



the other hand, the intervals must be wide enough to cover a 
large enough number of spectra. After some trials, we found that 
Az ~ 0.1 is a good choice. Hence, we produced SOMs for the 
quasar spectra in redshift bins with a 0.1 redshift step size. The 
mean number of spectra per bin amounts to 2 680, ranging from 
165 to 6 847. 



3.2. SOM algorithm 

The set of input data (i.e. spectra) is defined as vectors x(j) = 
[£i(j), ..,£nU)] e 2V 1 > where n = 975 is the number of pixels in 
each spectrum and j denotes the index in the sequence of source 
spectra j = . . . j max (where j max strongly varies with z). The 
neural network consists of i e {1 . . . N] neurons, represented by 
weight vectors mt = [//,- f i, . . . ,///,«] € % n , that are organised on 
a two-dimensional grid. We use a flat grid with closed bound- 
aries (in contrast to cylindrical or toroidal boundary conditions). 
The weight vectors are modified in each iteration step and thus 
vary with the discrete time coordinate t = 0, 1, 2, ... , £ max , i.e., 
nti = ntiit). Each weight vector ntt can be considered as an artifi- 
cial spectrum. The entire set of weight vectors approximates the 
distribution of input spectra. 

The SOM algorithm is essentially based on two processes 
that are responsible for the self-organising properties of the neu- 
ral network: (1) determining the best-matching unit (BMU) for 
a given randomly chosen input element and (2) the successive 
adaptation of the weight vectors in the neighbourhood of the 
BMU towards the given input element. 

For a given spectrum x(j), the Euclidean distance 
\\x(J) - nti\\ to each neuron m\ is computed and the BMU is iden- 
tified with the winning neuron m c , i.e. the neuron for which the 
distance is minimised, i.e., 



In the next iteration step (t — > t + 1), an adaptation of all neurons 
is performed, according to their difference from a spectrum Xj 
and the neighbourhood function h c j(r) 



niiit + 1) = nti(t) + h c4 (r) • [x(j) - mi(t)\ 
where r = t/t max . The neighbourhood function 



h cJ (r) = air) • exp 



/ lkc-r f || \ 
\ 2cr 2 (r) I 



(2) 



(3) 



acts as a smoothing kernel over the network and converges to 
zero with an increasing number of learning steps. The vectors 
r c e % 2 and r t e % 2 are the location vector of the BMU and 
and the location vector of the weight vector ntf, respectively. 
Compared to the frequently used Gaussian smoothing kernel, 
Eq. ^ uses a slightly modified version that has broader wings 
and a sharper peak. We found from various trials that Eq. ^ 
yields more reliable clustering results. 

The process of finding the BMU and adapting the weight 
vectors is performed for each input spectrum in random order 
within each learning step. The neighbourhood function is mod- 
ified over time by the learning rate a(r) and the radius func- 
tion cr(t) in such a way that (a) the map develops large-scale 
structures in the early phase, while (b) finer adjustments occur 
in later steps. Both functions are assumed to monotonically de- 
crease with time r € {0 ... 1} and are parametrised in the simple 
way 



a(r) = a I — J , cr(r) = cr \ — 
Wo/ Wo 



(4) 



IK/) - wdl = min,- {\\x(j) - will 



(1) 



where ao,ai,o~o, cr\ are the learning parameters of the Kohonen 
network with aro > &\ and cr > (T\ . The values used here are 
listed in Tab.Q] 

The number of neurons in the network must at least corre- 
spond to the number of source spectra but can be larger, i.e., 
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Table 1. Network parameters used for the clustering. 



Parameter 


Values 


Number of neurons N 


196... 8281 


Number of iteration steps t max 


100 


Learn radius cr 


1.0 


Learn radius cr\ 


0.125 


Learning rate a 


0.25 


Learning rate a\ 


0.01 



N > jmax- Our experience has shown that better results can be 
achieved when a certain fraction of positions within the map 
is not occupied with source spectra. These unoccupied neurons 
open up space to form cluster boundaries between distinct spec- 
tral types, which settle into clearly separated areas of the map. 
Outlier spectra have enough space to roam the neural landscape. 
We found that a ratio of N/j max « 1.2 represents a good trade- 
off. Hence, the grid size varies between 14 x 14 and 91 x 91 
neurons for the highest-z bin and the 1.5 < z < 1.6 bin, respec- 
tively. 

In the beginning of the iteration process at t = 0, each weight 
vector W/(0) is initialised with a random input spectrum x(j). 
The initialisation with purely random weight vectors would re- 
quire much more learning steps to obtain comparable results. 
In each iteration step, we start by computing the Euclidean dis- 
tances for the first input spectrum to all weight vectors within the 
network. The weight vector with the shortest distance is identi- 
fied with the BMU for this particular input spectrum. Hence, 
the input spectrum is moved to its network location for this par- 
ticular learning step. The BMU and all other weight vectors are 
then updated according to the neighbourhood function. This pro- 
cess is repeated for all remaining input spectra and over multiple 
learning steps until convergence is reached, i.e. until the BMUs 
and the associated input spectra essentially remain at the same 
positions in successive iteration steps. In our case of relatively 
small maps, the re-distribution of the spectra reaches sufficient 
convergence after about 100 iteration steps [j 

4. Selection and classification of unusual spectra 

4.1. Selection of spectra from the SOMs 

Originally, our search for unusual spectra was focused on un- 
usual FeLoBAL quasars. The longest- wavelength lines of the 
strongest UV iron multiplet, UV 1, (A < 2600 A) is shifted into 
the spectral window of the SDSS spectra for z ^ 0.6. This de- 
fines the lower limit for our quasar selection. The upper limit is 
set at z ^ 4.3 by the demands that (a) the Civ line should be 
well within the spectral window and (b) the number of quasars 
per z bin must not be too small. Therewith, the quasar sample is 
restricted to the range z = 0.6 . . . 4.3. Altogether, 99 151 spectra 
in this z range were analysed with the Kohonen method. 

For each of the 37 resulting SOMs, the corresponding grid 
of spectra was produced in the form of an "icon map", where 
the spectral resolution was reduced by an order of magnitude. 
Outlier spectra are clustered at the edges and corners of the 
SOMs (Sect.|3] Figs.[T] and [2]) and are thus easy to select by 
means of the visual inspection of the icon maps. The selected 



2 There are two reasons for such a small number of iterations in a 
neural network. First, our interest is focused on broadband spectral fea- 
tures and the spectral resolution is correspondingly low. Secondly, the 
final spectra maps are inspected visually. For an automated classifica- 
tion and higher resolution spectra, t max must be considerably increased. 



quasar sample is expected to contain various spectral types. We 
note, however, that the selection procedure (a) is certainly not 
free of a subjective bias, (b) does not use a quantitative crite- 
rion, and (c) is thus incomplete. The latter is, of course, a direct 
consequence of the term "unusual" not being defined a priori. 

We selected 1575 spectra of 1530 objects, 41 objects from 
this sample have more than one spectrum in DR7. The selected 
spectra were first corrected for Galactic foreground extinction 
and then individually checked to estimate z and object type. The 
Milky Way extinction curve from Pei (1992) was used for the 
extinction correction with E(B - V) obtained from the "Galactic 
Dust Extinction Service')^] of the NASA/IPAC Infrared Science 
Archive. 

While obtaining the redshift from a typical quasar spectrum 
is mostly straightforward, the situation is different for highly un- 
usual spectra. If prominent emission lines are present and clearly 
identified, the redshift is estimated as usual by comparing the 
measured line centres with a catalogue of their rest- frame wave- 
lengths. The spectral lines were selected manually to account for 
such effects as noise, broad absorption lines, or artefacts from 
bright night-sky lines. If the spectrum is not dominated by BAL 
troughs and emission lines are only rudimentarily present, or not 
at all, we tried to estimate the redshift by fitting the continuum 
to that of the SDSS quasar composite spectrum (VandenBerk et 
al. 2001 ). Redshift estimation was approached with particuliar 
caution for the BAL quasars without clearly identified emission 
lines but with many narrow absorption troughs or overlapping 
troughs. Here, we basically follow the approach outlined by Hall 
et al. (12002b . In several cases, redshifts could be estimated only 
roughly. Dubious as well as uncertain redshifts were flagged: 
f z = l for likely but not certain, 2 for uncertain due to the lack 
of redshift indicators, and 3 for uncertain due to low S/N, re- 
spectively. For example, z could not be estimated for three ob- 
jects with featureless blue spectra, which were flagged therefore 
fz = 2. 

4.2. Contamination 

The procedure of individually investigating each spectrum was 
repeated in three consecutive, independent runs. At the end, 
1014 of the selected 1530 objects (66%) were included in the 
(preliminary) catalogue of unusual quasars, 8% turned out to be 
stars, another 8% were classified as galaxies, and 5% were re- 
jected because of too low S/N and/or errors. About 13% were 
rejected as (nearly) normal quasars. Some of these quasars were 
positioned on the SOM in the area of unusual spectra because 
the assigned redshift from the SDSS pipeline was wrong. Others 
were simply located in transition regions between unusual and 
usual spectra but turned out to be, more or less, usual when in- 
dividually inspected. Objects displaying strong stellar continua 
and typical low-z absorption and/or narrow emission lines with- 
out broad components were simply classified as galaxies and 
subsequently rejected. Among the remaining 1014 catalogue en- 
tries, there are 18 objects for which z is very uncertain or not 
determined at all (f z > 1). 

Most of the rejected stars were of late spectral type, a few 
others are early-type stars and a small fraction were identified 
with exotic types. While normal stars were easily recognised 
based on their typical spectral features and rejected from the 
database, there remains a risk that some (a few) rather unusual 
stellar spectra were not correctly identified, e.g., rare types of 
white dwarfs (WD; Schmidt et al. 1999 2007; Szkody et al. 



http://irsa.ipac.caltech.edu/applications/DUST/ 
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Fig. 3. Examples for unusual spectra which were confirmed (panels a to g) or not confirmed (panel h), respectively, as rare types of 
stellar contaminants in the initially selected quasar sample. 



120041 Carollo et al. 120061) . Fig.[3]a-g show peculiar spectra that 
were classified as quasars by SDSS but rejected later on in the 
present study either during the inspection of the spectra or as 
a result of proper motion determinations (see below): (a) the 
superimposition of the spectra of a normal quasar and a late- 
type star (SDSS J014349.15+002128.3), (b) the WD- M star 
binary SDSS J100658.64+121133.9 (Heller et al. 120091 (c) a 
member of the new class of magnetic WD binaries with ex- 
tremely low mass tr ansfer rates (SDSS J103 100.55+202832.1; 
Schmidt et al. 2007), (d) a rare magnetic cataclysmic variable 
with extreme cyclotr on features (SDSS J13241 1.57+032050.4; 
Szkody et al. 2004), (e) the newly discovered proper motion 
star SDSS J02073 1.81 +004941. ^(Table^. In a first attempt, 
we had classified SDSS J02073 1.8 1+00494 1.6 as a low-redshift 
FeLoBAL quasar (similar to SDSS Jl 12526.12+002901.3 and 
SDSS Jl 12828.31+01 1337.9 from Hall et al. 2002). Its signif- 
icant proper motion suggests that SDSS J02073 1.8 1+00494 1.6 
is a possible extreme DZ or carbon star, (f) the newly discov- 
ered proper motion DQ WD SDSS J020534.13+215559.7, and 
(g) the featureless blue spectrum of a newly discovered proper 
motion star (SDSS J084357.71+115624.1). 

Finally, SDSS J014733.58+000323.2 (Fig.[3]h) is an exam- 
ple of a spectrum that was erroneously rejected in our first anal- 
ysis attempt. The object was classified as a stellar mass black- 
hole candidate by Chisholm et al. (2003). Though there is a blue 
continuum with apparent broad emission-line components, these 
lines do not fit the positions of the typical broad lines in quasar 
spectra. Our referee, Dr. Patrick B. Hall, pointed out that there is 
Civ/Li 1548.20,1550.77 absorption at z=1.8815. This high red- 
shift clearly proves that SDSS J014733.58+000323.2 is a quasar. 
The closer inspection of the spectrum leads to its classification as 
a slightly reddened (E(B - V) = 0.07 for SMC-like dust) weak- 
line quasar at z = 1.98 without Cm] emission and with a shift 
between the C iv and Mg n emission lines. 

Stellar contamination remains a serious problem even after 
the individual inspection of all spectra. Various studies have used 
astrometric information for the quasar selection (e.g., Sandage & 
Luvten lT967l Kron & Chiu flMTt Meusinger et al. 12002120031 



Kaczmarczik et al. 2009; Lang et al. 2009). Quasars are suf- 
ficiently distant to able to neglect there having any measur- 
able absolute proper motion (pm) in presently available posi- 
tion data. Therefore, absolute pm are considered helpful data to 
discriminate quasars from (nearby) stars. We cross-correlate our 
quasar sample with the PPMXL catalogue (Roser et al. 2010). 
This catalogue is the largest collection (of nearly one billion 
objects) of pm in the International Celestial Reference System 
(ICRS), which is primarily realized by the Hipparcos catalogue. 
The PPMXL aims to be complete down to V ~ 20, while our 
quasar sample extends to fainter magnitudes (mean g magnitude 
(g) = 20.02 + 1.41). With a search radius of 5", we identified 
918 objects (91%) from our sample in the PPMXL (mean po- 
sition difference 073). We express the probability of significant 
non-zero pm in terms of a simple pm index, which is defined as 
the total pm \i in units of the pm error e\ 



/ P m = ^(Ma cos S) 2 +ti 2 6 j ^Je 2 aCOs6 



+ e, 



V6 ' 



(5) 



For the 918 identified objects, we found a mean value of (7 pm ) = 
2.46 and a median at 0.98; 90% of the objects have 7 pm < 2.5. We 
restricted the further analysis to the 95 objects with 7 pm > 2.5. 
One third of this subsample have large pm [i > 150 arcsec yr" 1 , 
which are likely due to these objects beeing fakes (Roser et al. 
2010). We re-investigated the spectra of all 95 pm candidates 
and searched for possible entries in the SIMBAEJj and NEEJj 
database. For only one object, SDSS J020534.13+215559.7, did 
the spectrum turn out to be characteristic of a DQ WD (Fig. [3]). 
The star was not found in any stellar database. For 81 objects, 
the spectra could be assigned to quasars without doubt, another 
9 objects are likely quasars. Another four selected pm candi- 
dates have relatively featureless spectra, one of them, SDSS 
J 153939. 1+274438, is a known radio source with dominant 
blazar characteristics (Massaro et al. 2009 ). The large AGN frac- 
tion among the pm- selected objects obviously means that the 
accuracy of the available pm data is insufficient to efficiently 
select stellar contaminants in the magnitude range of our sam- 
ple. The trend of increasing 7 pm with fainter magnitudes (e.g., 



4 According to the SDSS explorer, a quasar with z 
and z = 2.3273 in DR8, respectively. 



1.1165 in DR7 



http://simbad.u-strasbg.fr/simbad/ 
http://ned.ipac.caltech.edu 
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</ pm ) = 1.12 for g < 20 compared to </ pm > = 4.10 for g > 20) 
and vice versa (e.g., (g) = 19.7 for 7 pm < 2.5 compared to 
(g) = 21.0 for 7p m > 2.5) indicates that the pm errors are sys- 
tematically underestimated at faint magnitudes. 

In a second step, we therefore used all available multi-epoch 
positions from different sky surveys to determine an improved 
pm. The longest time baseline (with epochs between 1950 and 
2002) was provided by the SuperCOSMOS Sky Surveys (SSS; 
Hambly et al. 2001 ), including the measurements of overlap- 
ping Schmidt plates. In the case of problems with the SSS mea- 
surements on the first epoch Palomar Observatory Sky Survey 
(POSS1) plates, we used the APM measurements of these 
POSS1 plates (McMahon et al. 12000ft . New data (for epochs be- 
tween 1998 and 2009) came from two data releases of the SDSS 
(DR7, Abazajian et al. 2009; DR8, Aihara et al. 2011; with 
30-70 epochs in the equatorial stripe), the eighth data release 
of UKIDSSQ and for the brightest objects, from the Carlsberg 
Meridian Catalog (2006), and the Two Micron All Sky Survey 
(2MASS; Skrutskie et al. 12006ft . Finally, we included the epoch 
2010 positions of the objects (if measured) in the WISE (Wide- 
field Infrared Survey Explorer; Wright et al. 12010ft preliminary 
data release. The pm obtained from simple linear fitting of all po- 
sitions over time were not only judged in terms of their formal 
errors (typically smaller than the PPMXL errors, if available), 
but also by the influence of the less-accurate Schmidt plate data. 
In particular, we did not trust any formally significant pm if the 
fit was dominated by only one POSS1 position and could not be 
confirmed after excluding this plate from the fit. Examples of the 
proper motion fits are shown in Fig. [4] 



SPSS J020534+215600 




SDSS J151627. 4+30S220 



1980 
Epoch 



Fig. 4. Two examples for estimating proper motions of uncertain 
quasar candidates: the DQ WD SDSS J020534. 13+215559.7 
(top) and the quasar SDSS J151627.40+305219.7 (bottom). 

The procedure was applied to objects from three subsamples: 
(a) the 95 pm candidates from the PPMXL, (b) the two "myste- 
rious" object s fro m Hall et al. (2002) and 15 possibly related 
objects (Sect. |6.2| ), and (c) an additional 9 objects that had ei- 
ther featureless spectra or spectral features we were unable to 
interpret. Among these 121 objects are in particular all those 
with f z = 2. As expected, most quasar candidates have zero 
proper motions that we were able to measure very accurately if 
many epochs were available. However, for seven candidates we 
found significant proper motions (Table[2]). With the exception 
of SDSS J144354.73+242906.6, the spectra of all these objects 



are blue and suggest that their classification as WDs with the 
spectral types given in the last column of Tab.|2] The spectra of 
SDSS J134913.51+205646.9 and SDSS J144354.73+242906.6 
are noisy; the former is probably a DC and the latter may be a 
late-type star, but the S/N is too low for a reliable type estimation 
to be possible. Another two unusual quasar candidates, SDSS 
J032907.24+002438.4 and SDSS J033716.08+000041.7, turned 
out to be galaxies. These 9 entries were removed so that our 
catalogue finally contains 1005 entries, among them 14 objects 
with f z >2 that were excluded from the statistical investigations 
described in Sect.[5j 

Table 2. Rejected quasar candidates with significant proper mo- 
tions from the present study (N Q \ number of epochs, TVd: WD 
spectral type). 



SDSS J 


\i a cos 6 
(mas/yr) 


(mas/yr) 


N e 


Twd 


020534.13 + 215559.7 


109.0 + 7.2 


+ 10.3 + 5.6 


5 





020731.81+004941.6 


29.2 + 1.0 


-6.2+1.5 


39 


Z 


024058.80 - 003934.5 


-10.4 + 0.9 


-7.5 + 1.2 


80 


C 


084357.71 + 115624.1 


-9.4 + 8.0 


+45.8 + 6.0 


9 


C 


100149.22 + 144123.8 


-345.5 + 1.5 


-2.3 + 3.5 


7 


C 


134913.51+205646.9 


-68.4 + 8.7 


-17.6 + 3.8 


5 


? 


144354.73 + 242906.6 


-107.7 + 6.2 


-112.8 + 7.1 


6 


? 



The UKIDSS project is defined in Lawrence et al. (2007 ) 



As emphasised above, our method of selecting unusual 
quasars is not aimed at a complete sample. Nevertheless, to 
get an idea of the completeness we considered the sample of 
the unusual BAL quasars from Hall et al. (2002). Among these 
23 objects, four were classified by the SDSS DR7 as unknown 
(spec_cln=0), and for another three the spectra were unavailable 
for download. Hence, 16 quasars from this sample are expected 
to be selected by our approach. Only one Hall quasar is missed, 
SDSS J121441.42 - 000137.8, a BAL quasar with possible rela- 
tively strong Fe in absorption. 

4.3. Categories of unusual quasar spectra 

According to the different areas populated by unusual quasar 
spectra in the Kohonen maps, a classification into essentially 
four categories was immediately suggested. The corresponding 
dominant spectral features are (a) strong BALs, (b) red continua, 
(c) weak or absent emission lines, and (d) strong optical or UV 
iron emission. As a result of evaluating all selected spectra, we 
decided to subdivide the first category into two types: unusual 
BAL structures (mostly LoBALs) and more or less normal BAL 
quasars (mostly strong HiBALs). Moreover, we found that the 
category of red quasars contains a substantial fraction of spec- 
tra where the continuum is significantly red or reddened at UV 
wavelengths (e.g., A ^ 3000 A), but seem to reasonably fit the 
SDSS quasar composite spectrum at longer wavelengths (see 
also Hall et al. 2002] their section 6.3.1). Therefore, we split 
this category into the two types of (pure) red continua and UV- 
red continua, though it was unclear ab initio whether they rep- 
resent distinct types. Finally, a relatively small number of ob- 
jects could not be adequately sorted into one of these six types. 
These quasars form a small and quite inhomogeneous group, in- 
cluding those that have spectra with exceptionally blue contin- 
uum, unusual line profiles, unusual continuum shape (perhaps 
caused by an error in the flux calibration), or quasars with nar- 
row emission lines but lacking substantial broad emission-line 
components (type 2 quasar candidates). 
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Fig. 5. Example spectra for the unusual quasar types A to M. For comparison, the SDSS quasar composite spectrum (VandenBerk 
1200 II) is shown, arbitrarily normalized to each object spectrum at the red end. 



Therewith, this paper deals with seven categories or types of 
unusual quasar spectra (Fig. [5}: 
A - Unusual BAL quasars (mostly LoBAL quasars); 
B - Quasars with strong but more or less usual BAL structures 

(mostly HiBALs); 
C - Quasars where the whole continuum is significantly redder 

than the SDSS composite; 
D - As type C, but strong reddening is obvious only in the UV; 
E - Quasars with weak or absent emission lines; 
F - Quasars with strong iron emission; 
M - Miscellaneous. 

Since some of the selected spectra appear to be very complex, we 
decided to build the classification on the subjective impression 
of the overall picture rather than on quantitative criteria. 

The transitions between the various categories are smooth. 
Many spectra show a mixing of characteristic features of two 
or three types (see e.g., Fig.[?]). Therefore we allocated to each 
quasar up to three out of the seven possible types T- u i = 1 ... 3 
(T\ ± 72 ^ T$) ), where the significance decreased with i. Forty- 
one percent of the quasars have more than one classification and 
8% have more than two. In particular, there is a strong coupling 
between unusual BALs, red continua, and strong iron emission. 
Table[3] lists all three types. In the rest of the paper, we consider 
the type of highest significance (T = T\) only. 

4.4. The catalogue 

The results of our visual inspection of the spectra of the 1005 
unusual SDSS quasars are given in Table[3] The catalogue was 
matched with both the QCDR7 (Schneider et al. 120101) and the 
08Jull6 version of the catalogue from the FIRST Survey (Becker 
et al. [T995). The vast majority (98.3%) of the quasars have en- 
tries in the QCDR7 (quasar catalogue flag f c > 0). The agree- 
ment between the redshifts from the present study and those 
from the QCDR7 is generally goodjjAs expected, the discrep- 
ancies are larger when the redshifts from the SDSS pipeline are 

8 For the quasars with certain redshifts (f z = 0), we found six with 
deviations > 10%. These objects are indicated by f q = 2. The devia- 
tions were larger than 50% for SDSS J093437.53+262232.6 and SDSS 
J131524.00+041734.4. 



considered: 98 (41) quasars with deviations > 10(50)%. The 18 
objects missed in the QCDR7 have normal absolute magnitudes 
and redshifts (M{ = -25.5, z = 1.84) but their redshifts are 
mostly uncertain (6/7/4/0 objects have f z = 0/1/2/3). For these 
objects, the entries in cols. 6 and 8 of the catalogue were taken 
from the SDSS DR7 Explorer. 

The full catalogue is available only in electronic form. 
Table[3] lists the first five entries for guidance regarding its 
content and form: 

Column (1): The running catalogue number. 
Column (2): The SDSS J2000 equatorial coordinates for the 
quasar (taken from the spectrum fits header). 
Column (3): The QCDR7 catalogue flag (f c = 0: not in QCDR7, 
1: in QCDR7 with z deviation < 10%, 2: in QCDR7 with z 
deviation > 10%). 

Columns (4) and (5): The redshift z and the redshift flag f z 
from the present study (f z = 0: certain, 1 : likely but not certain, 
2: lack of certainly identifiable redshift indicators, 3: low 
signal-to-noise). 

Column (6): The SDSS apparent i band magnitude corrected for 
Galactic foreground extinction. 

Column (7): The absolute i band magnitude, computed as 
described by Kennefick & Bursick (2008 ) for a spectral index 
a v = -0.5. 

Columns (8): The FIRST radio detection flag from the QCDR7 
(0 for non-detection). 

Column (9): The FIRST 1.4 GHz peak flux F lA (mJy). 
Column (10): The radio loudness parameter R{ computed from 
F1.4 and the extinction-corrected SDSS i band magnitude 
following Ivezic et al. (120021) . 

Column (11): The peculiarity index x 2 °f the spectrum 
(Sect.|5J). 

Columns (1 2) to (14): The type classifications from the present 
study (Sect. |4.3| ) in numerical format (1=A, 2=B, ...) for com- 
putational reasons. The priority decreases from T\ to T$. 
Column (15): Comments on spectral peculiarities. 
Abbreviations: 

- abs: absorption 

- ao: associated absorption 

- bel: blueshifted emission line 



10 



H. Meusinger et al.: Kohonen- selected unusual SDSS quasars 



Table 3. Classification and properties of unusual SDSS quasars. (Only the first five rows are shown here; the full table is available 
electronically only.) 



Nr 


SDSS J 


/c 


z 


fz 


k 


Mi 


/r 


Fia 


Ri 


x 1 


T\ 


T 2 


T 3 


Comment 


1 


000009.38 + 135618.4 




2.233 


1 


18.17 


-27.45 









2 


5 


6 




Fe in em? 


2 


000326.67 + 001157.1 




1.845 





18.95 


-26.23 









18 


4 








3 


000625.25 + 154625.8 




1.527 


1 


19.58 


-25.16 









77 


3 






ao(Mg n); strong [Ne v];noisy 


4 


000728.45 - 042345.5 




1.198 





18.12 


-26.04 









15 


6 


1 




Fein em; see 2215(H02) 


5 


000920.01+005618.4 




2.387 





20.20 


-25.56 









10 


4 


5 







- bq: blue quasar 

- em: emission 

- dlas: damped Ly a absorption system 

- fg: foreground 

- 1st: longwards-of-systemic trough 

- mnt: many narrow troughs 

- my st?: related or possibly related to the two mysterious ob- 
jects from Hall et al. (l2QD2l) 

- nlq: narrow-line quasar 

- nt: narrow troughs 

- ot: overlapping troughs 

- sdvdpc: spatially distinct velocity-dependent partial covering 
of the continuum source 

- see 0810(H02): see object SDSS J081024.75+480615.5 
from Hall et al. (12002b 

- uc: unusual continuum 

- ulp: unusual line profile 

5. Mean properties of the various unusual quasar 
types 

The number of quasars in the seven types are listed in Table|4] 
along with some mean properties discussed below in this section. 



Table 4. Mean properties of the unusual quasar types. 



T 


N 


x 2 


z 


M { 


/-(intr) 


M (cor) 

i 


/rl 


A 


215 


78.8 


2.02 


-26.3 


0.18 


-27.6 


0.26 


B 


209 


12.7 


2.92 


-27.0 


0.09 


-27.9 


0.11 


C 


116 


61.3 


1.96 


-25.9 


0.38 


-28.5 


0.31 


D 


150 


115.6 


1.47 


-25.9 


0.20 


-27.0 


0.30 


E 


185 


9.3 


2.14 


-27.0 


0.00 


-27.0 


0.26 


F 


112 


28.5 


1.40 


-26.2 


0.06 


-26.5 


0.16 


M 


18 


13.0 


1.53 


-25.7 






0.33 


all 


1005 


48.8 


2.06 


-26.4 


0.15 a 


-27.4 a 


0.23 



Notes. (a) Types A to F only. 



5.1. Peculiarity index 



As a quantitative measure of the deviation of a given spectrum 
from the SDSS quasar composite spectrum ( VandenB erk 1200 1 1) . 
we computed 



X 



1 N 

= /V^2 Zj [ FnWi) " F comp,n(^i)] 



(6) 



i=l 



where F n and F com p jn are the normalised spectra of the quasar 
and the composite, respectively, and cr is the average noise of 



F n . The wavelengths A{ refer to the observer frame and the index 
i indicates the pixel number. The quasar spectrum is normalised 
to the integrated flux with an integration interval from 4000A 
to 9000A for z < 2.29 and from 1216 • (1 + z)A to 9000 A for 
z > 2.29, respectively. The SDSS quasar composite spectrum 
was shifted to the redshift z and normalised to match the red end 
of F n (A) dXX- 8600-9 100A. The noise cr was derived from the 
difference spectrum D n (A) = F n (A)-F n ^io(A), where F n , s io(^) is 
the smoothed version of F n (A) using a 20 pixel boxcar filter. To 
reduce the effect of residuals from narrow spectral features, five 
rest-frame wavelength intervals of ~ 50A width were selected 
where the contribution from emission lines is small (pseudo- 
continuum windows; see Sect. |5.5| ). The average standard de- 
vation of D n (A) in these parts of the spectrum measured in the 
observer frame is taken as a proxy for cr. Care has been taken to 
exclude the regions around the strongest telluric emission lines 
at ^5577A and /16300A. 

At least for the statistical comparison of subsamples, x 2 can 
be taken as a useful integral measure of the "peculiarity" of the 
spectrum. The individual values cover the broad range from ~ 1 
to 5 000 with a mean value of ~ 50. For comparison, the 15 un- 
usual Hall quasars (Hall et al. 120021) cover the range from 35 to 
1089 with a mean value of 261. When we restrict the selection 
to x 2 > 35, we have 296 quasars corresponding to an unusual 
quasar fraction of / uq = 296/99151 = 0.0029 in perfect agree- 
ment with 23/8000 = 0.0029 for the selection by Hall et al. 
However, if we choose the threshold so that our selected sample 
has the same mean value of x 2 as the Hall sample, our selection 
is much less efficient with / uq = 0.0012. 

The peculiarity index x 2 is particularly sensitive to devia- 
tions from the composite spectrum over wide wavelength inter- 
vals, i. e., owing to intrinsic reddening or an intrinsically red 
continuum. It is thus unsurprising that the largest values are 
measured for the red/UV-red quasars of the types C and D, but 
also type A. (However, type C spectra frequently suffer from 
strong noise, i.e., large values for cr). For weak-line quasars, on 
the other hand, the deviations are restricted mainly to the small 
wavelength intervals of the lines and thus x 2 is generally small. 
The unusual BAL quasars of type A have significantly larger 
values of^ 2 than their more normal relatives of type B, in agree- 
ment with their stronger apparent intrinsic reddening (see be- 
low). 

5.2. Redshifts and absolute magnitudes M- x 

The first row of Fig. [6] shows the histograms of the z distribu- 
tions. As can be seen, the different types cover different redshift 
intervals. This is simply due to the limited spectral window of 
the observations in combination with the characteristic spectral 
features being tied to special wavelength intervals. For example, 
the change in the continuum slope for the type D quasars fre- 
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Fig. 6. To/?: Redshift distribution for types A to F. Second and third row: absolute i magnitude and 1.4 GHz luminosity, respectively, 
as function of redshift compared with QCDR7. Bottom: Relative distribution of the 1.4 GHz luminosity. 



quently appears around rest- frame wavelengths between 2000A 
to 3000 A and is thus best observed f or z ~ 1.5. At higher z, the 
wavelength of the turnover is shifted towards and beyond the red 
edge of the spectral window and the spectrum appears to be red 
(see also Fig. [5]). Strong Fe emission in the UV is restricted to 
the same wavelength interval. Consequently, quasars of types D 
and F cover similar z intervals. 

In the second row of Fig.[6j the absolute magnitudes M{ are 
plotted as a function of z (crosses), along with the median re- 
lation and l<x deviations for the QCDR7 quasars. Owing to the 
strong inherent redshift bias, the M{ distribution of any type can- 
not be compared with that from the whole QCDR7. For the same 
reason, the M- x distributions of the various quasar types cannot be 
compared with each other. Therefore, we constructed for each 
type corresponding comparison samples from the QCDR7 that 
have identical z distributions. For every quasar (redshift z, flag 
f z < 2), we randomly selected one quasar from the QCDR7 with 
a redshift in the interval z - 0.05 . . . z + 0.05. This procedure was 
performed 100 times to create 100 different comparison samples. 

To test whether our unusual quasars and their comparison 
samples represent the same quasar population with respect to 
Mi, we applied the two-tailed two-sample Kolmogorov-Smirnov 
(KS) test (e.g., Siegel & Castellan [W88]) . The null hypothesis 
HO, that the members of a given type of unusual quasars have 
the same M- x distribution as the quasars in a comparison sample 
from QCDR7, is tested against the alternative HI, that the M- x 
distributions of the two samples are different. The KS test uses 
the maximum difference D = max/ \S t - S^ comp \ between the 
cumulative distributions S{ from the two samples. For the two- 
tailed test, HO has to be rejected at a chosen level a if D is so 



large that the probability p of its occurrence is p < a. Since we 
have 100 comparison samples, we performed 100 tests per type. 
Adopting a = 0.05, HO had to be rejected for all comparison 
samples of types E and F. For the other types the null hypothesis 
has to be rejected for only 45% of the tests for type C and for 
< 10% for the types A, B, and D. 

According to Fig. [6] the quasars of types E and F tend to 
be more luminous than normal. We applied a one- tailed KS 
test with the alternative hypothesis HI that unusual quasars are 
more luminous than usual quasars. Again, H0(a = 0.05) had 
to be rejected in favour of HI for the types E and F. We con- 
clude that both the weak-line quasars and the strong iron emit- 
ters are more luminous in the optical than typical quasars from 
QCDR7. Just et al. (2007) created a sample of 32 of the most 
luminous quasars from the SDSS DR3 quasar catalogue adopt- 
ing a minimum luminosity at M- x - -29.28. Our whole sam- 
ple includes three quasars brighter than this threshold; all three 
are of type E, among them the ultra-luminous quasar SDSS 
J152156.5+520238 with M { = -30.17. 

It is conceivable that our selection process for spectra with 
weak emission lines or with strong iron emission might have 
preferentially selected quasars with higher spectral S/N, which 
will be the more luminous ones at each redshift. We estimated 
S/N for the unusual quasar spectra as well as for the spectra of 
~ 10 4 comparison quasars. The noise was measured in defined 
areas of the spectrum as described in Sect. |5.1| the signal was 
given by the mean flux in the same areas. We indeed found a 
higher average S/N for the unusual quasars namely 11.3 and 15.6 
for types E and F, respectively, compared to 8.1 and 9.2 for the 
corresponding comparison samples. This is, however, not com- 
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Fig. 7. Colour indices (extinction corrected) as function of redshift for types A to F (top to bottom) in comparison with the SDSS 
DR7 quasar catalogue (single values and median relations). 



pletely unexpected because of the trend of S/N with apparent 
magnitude. In the next step, we created new comparison sam- 
ples where for each unusual quasar one comparison quasar was 
identified with similar z (as before) and with similar S/N. For 
either type E and F, we found that the absolute magnitudes M- x 
of the unusual quasars of types E and F are on average about 0.4 
mag brighter than the comparison samples of normal quasars. 
The KS test confirms that the differences are significant. 

Among all six types, the red quasars of type C have the 
lowest optical luminosities, as expected. We performed the one- 
tailed KS test to check whether the red quasars are significantly 
underluminous and found that this is unlikely to be the case. 

5.3. Radio luminosity and loudness 

The radio detection fraction, / RD , of our whole sample is 0.29, 
compared to 0.08 for the quasars from the QCDR7 in the same 
redshift range (0.5 < z < 4.4). A similar difference is seen for 
the fraction /rl of radio-loud quasar^j which is 0.22 in our sam- 
ple, compared to 0.07 for the QCDR7 quasars. This could imply 
that independently of colour selection, unusual quasars are de- 
tectable over a much longer time span at the FIRST level than 
normal quasars. Alternatively, these significant differences may 
be caused by the SDSS having also targeted FIRST sources for 
spectroscopy when their optical colours did not meet the quasar 
selection criteria. The FIRST target flag is indeed set for 21% of 
the quasars in our sample, compared to 5% in the whole quasar 
catalogue. 



9 Using the criterion R[ 
Ivezic et al. ([20021 . 



> 1 as definition for radio loudness; see 



To test whether the excess radio-detected unusual quasars 
were only discovered through FIRST targeting, we considered 
the numbers N? of quasars with the FIRST target flag set and Nq 
of quasars with their colour target flag set. We found N^/Nq = 
0.33 for the unusual quasars compared to 0.07 for the QCDR7. 
A similar result was found when we compared the number Afps 
of quasars selected solely by the FIRST selection (but not the 
colour selection) to the number A^c s of quasars without FIRST 
counterparts but selected solely based on their colours, namely 
AWAfcs = 0.22 for the unusual quasars and 0.01 for the QCDR7. 
Moreover, when we binned the unusual quasars into intervals 
of x 1 , it could be clearly seen that the fraction of quasars with 
their FIRST target flag set increases with the mean peculiarity 
index (Tab.|5). Both N F /N C and N F JN Cs rise strongly with^ 2 
and reach values of ~ 1 for the objects showing the strongest 
mean deviation from the SDSS quasar composite spectrum. We 
could identify no physical reason for such a trend and concluded 
therefore that the increase in the fraction /rd of quasars with 
FIRST detections towards the more peculiar spectra is more 
likely a selection bias: a large fraction of the unusual quasars 
were not selected by the colour selection criterion but were tar- 
geted by SDSS just because they had been detected as FIRST 
radio sources. 

A similar trend, but weaker, was also observed for the frac- 
tion /rl of the radio-loud quasars. On the other hand, no indi- 
cation of this trend is seen for the ratio /rl//rd- The fraction 
of radio-loud among the radio-detected quasars for our whole 
sample is 0.76, compared to 0.90 for the QCDR7 sample. This 
means that our selection criterion of peculiar spectra did not in- 
duce a substantial bias towards quasars that are more active at 
radio frequencies than normal. 
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Table 5. Radio detection fraction for different spectral peculiar- 
ities. 



x 1 


N 


N F /N C 


N F JN Cs 


/rd 


/rl 


/rl/Trd 


0...5 


248 


0.17 


0.09 


0.19 


0.18 


0.94 


5 ... 20 


333 


0.20 


0.15 


0.23 


0.16 


0.72 


20. ..50 


184 


0.34 


0.22 


0.28 


0.23 


0.82 


50... 100 


123 


0.69 


0.41 


0.48 


0.38 


0.80 


100... 200 


61 


0.73 


0.62 


0.48 


0.39 


0.83 


>200 


42 


1.16 


0.91 


0.55 


0.31 


0.57 



the catalogue of BAL quasars from Trump et al. ([2006). Shankar 
et al. (2008 ) argue that such a trend fits well within a simple ge- 
ometric model. Together with the results in Tab.[5] the increase 
in / towards lower L\ A for type A indicates that the real fraction 
of unusual BAL quasars might be considerably larger than in 
present samples, where these quasars have instead been targeted 
mainly as FIRST radio sources. 

Finally, we note that there is a similar trend of /(log L\ A ) for 
types B and C. This trend seems to also exist for type F but is 
not statistically significant. 



On the other hand, it is interesting to study whether our lower 
radio loudness fraction is significant and how it depends on the 
type properties. For this purpose, we converted the FIRST radio 
flux densities F\ A to the specific luminosities L\ A emitted at 1.4 
GHz (restframe) via 

logL L4 = 23.08 + 2 log D L + log F 1A - (1 + a) log (1 + z\ (7) 

where L\ A and F\ A are given in WHz -1 and mJy, respectively, 
D L is the luminosity distance in Gpc, and a is the spectral index 
(F v oc y '), where we assumed that a - -0.5. The results are 
shown in Fig. [6] for the quasars of each type (crosses). The (nor- 
mal) quasars from the QCDR7 are represented by the median 
relation and the area populated by 80% of the quasars closest to 
the median on either side. 

As discussed above for the absolute magnitudes, comparison 
samples with the same z distribution are needed to check whether 
the radio luminosities in our sample differ from those of normal 
QCDR7 quasars. Such comparison samples were constructed in 
the same way as for M- u with the only exception that now only 
the subsamples of the radio-detected quasars (F\ A > 1 mJy) 
were considered. We again first performed the two-tailed KS test 
with the null hypothesis HO, that there is no difference between 
the radio luminosities of the unusual quasars of a given type and 
the comparison sample from the QCDR7, and the alternative HI, 
that both groups have different radio properties. We found that 
HO has to be rejected for all comparison samples of type A and 
for 80% of the comparison samples of type F. For the other types, 
the fraction of rejections is 20% or less. Since the fraction of 
radio-loud quasars tends to be smaller for the unusual quasars, 
we also applied one-tailed tests with HI : the radio luminosities 
of the unusual quasars are lower. We found that HO has to be 
rejected in favour of if 1 for the unusual BAL quasars of type A 
(p = 0.01) and the strong Fe-emitting quasars (p = 0.04) of type 
F. No such firm conclusion can be drawn for the other types. 

The bottom row of Fig. [6] shows the distribution of the ra- 
tio / = N(Li A )/Nqc(Lia) where N(Li A ) is the differential radio 
luminosity distribution for the corresponding type with equidis- 
tant logLi.4 intervals and Nqc(L\ A ) is the distribution averaged 
over 10 comparison samples from QCDR7, where both distri- 
butions are normalised. The vertical bars indicate the propaga- 
tion of the counting errors of N(L\ A ) and Nqc(L\ A ), which were 
identified with the Poisson limits at a confidence level of 0.84^1 
To compute the upper and lower Poisson levels, the approxi- 
mate equations given by Gehrels (1986) were applied. It can be 
clearly seen that strong BALs are not among the most radio-loud 
quasars. The fraction of BAL quasars drops by a factor of at least 
three between log L\ A = 25 and 27. This decline in the fraction 
of BAL quasars with increasing radio power confirms previous 
results from Becker et al. (2000) for a smaller sample of BAL 
quasars and Shankar et al. (2008 ) for a larger sample based on 



5.4. Optical colours 

Systematic trends of the SED can be illustrated by average spec- 
tra but also by the distribution of the colour indices. We first 
discuss the colours. 

Since the colour indices of quasars vary with z, it is use- 
ful to study colour-redshift diagrams rather than colour-colour 
diagrams. Such colour-z diagrams are displayed in Fig. [7] for 
u- g,g - r,r- i, and r - K, where u, g, r, / are the SDSS mag- 
nitudes from the QCDR7 corrected for Galactic foreground ex- 
tinction and K is the K band magnitude from the 2MASS survey 
(Skrutskie et al. 12006b . About 60% of our quasars have 2MASS 
measurements. We do not see significant differences between the 
subsamples with or without 2MASS counterparts. As in Fig.[6j 
the quasars from the QCDR7 are plotted for comparison (median 
and l<x deviation). 

In general, BAL quasars are known to be significantly red- 
der than those without BALs, and LoBAL quasars are even red- 
der (Sprayberry & Foltz [T9921 Reichard et al. [20031 Gibson 
et al. [2009). n can fo e clearly seen in Fig. [7] that types A to D 
tend to have redder colours than the whole SDSS quasar popula- 
tion with strongest deviations for type A. The colours of type C 
and D quasars are by definition redder than normal. For type E, 
the colour-redshift diagrams resemble those of the whole quasar 
population. The diagrams for type F quasars indicate a moderate 
reddening of u - g and g-r, while r - i and r - K are approxi- 
mately the same as for normal quasars. 



Corresponding to Gaussian statistics Icr 



As shown in Sect. |5.2| types A to D have the same absolute 
magnitudes as normal quasars. However, if their red colours are 
due to intrinsic reddening by dust and/or gas, the intrinsic lu- 
minosities of these quasars appear to be higher than normal. In 
Sect. |5.6| we estimate the intrinsic absorption and perform the 
KS test for the correspondingly corrected M- x . 

5.5. Composite spectra 

Composite spectra are particularly useful for determining the av- 
erage spectral properties of quasar samples (Francis et al. 1991 ; 
VandenBerk et al. 2001; Richards et al. [2003]) . The general 
procedure of constructing quasar composites is the following: 
after the redshifts are determined, the individual (foreground) 
extinction-corrected spectra are rebinned to a common wave- 
length scale in their restframe, and the quasars are sorted by 
redshift. Then, generating composites requires essentially two 
steps: (a) normalisation of the spectra and (b) combining the nor- 
malised spectra. Both steps can be performed in different ways. 
The normalisation is not entirely trivial because the spectra 
cover different rest- frame wavelength intervals. To create global 
composite spectra of normal SDSS quasars, Vanden Berk et al. 
(2001) started with the arbitrarily scaled spectrum of the lowest 
redshift quasar (with measured [O in] A5001 line emission), and 
scaled each subsequent spectrum by the overlap of the preceding 
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Fig. 8. Arithmetic median composite spectra for types A to F. For comparison the SDSS quasar composite spectrum from 
VandenBerk et al. (2001 ) is shown, arbitrarily normalized at the red end. The dashed vertical lines indicate the strongest absorption 
lines for types A and B. 



average spectrum (i.e., by its mean flux density in the wavelength 
interval of the overlap). The same procedure was applied in our 
previous study (Meusinger et al. 12011b . Here, however, we fol- 
lowed a slightly different approach. The main reason is that, for 
our present sample, the flux density of the overlapping interval 
can be strongly modified by the spectral peculiarities for which 
we search. Hence, we decided to restrict the normalisation to a 
few, relatively narrow continuum windows. Though there is no 
ideal continuum window where a suppression by BALs and/or 
the contribution from emission lines can be completely ignored, 
usable pseudo-windows exist however at 3540-3600A, 3030- 
3090A (Tsuzuki et al. 120061 . and 2000-2300A (Sameshima et 
al. 120111) . We started by scaling the lowest-redshift spectra. If 
spectra existed with z < 1.5, they were arbitrarily normalised 
to an average flux density F A = 1 over the wavelength interval 
A = 3540-3600A. The thus normalised spectra were combined 
to a first composite C\ . Next, the spectra in the redshift range 
z =1.5-1.94 were scaled to the flux density of C\ at A - 3030- 
3 090 A and the normalised spectra were co-added together with 
the normalized spectra from the previous step to create a second 
composite C^- In the same way, the continuum window at A - 
2200-2230A was used to extend the composite C3 to z =3. To in- 
clude quasars at higher redshift, we introduced another pseudo- 
continuum window at A - 1650- 1700 A where the spectra with 
Z > 3 were scaled to the flux density of C3 . 

Combining the normalised spectra to a composite requires 
choosing between mean or median and arithmetic or geometric 
averaging. Following Vanden Berk et al. (120011) . we used com- 
bining techniques to create (a) the arithmetic median, which pre- 
serves the relative fluxes along the spectra, and (b ) the geometric 
mean spectrum, which preserves the global shape of a power-law 
continuum. The arithmetic median spectrum was computed in 



the same way as in our previous study (Meusinger et al. 1201 II) : 
All normalised (restframe) spectra S t of a given type were in- 
serted into a 2D image, one spectrum per row with 1 A binning 
on the horizontal axis. The arithmetic median was then com- 
puted by averaging over each column using the procedure av- 
erage/image from the ESO-MIDAS packagq^] with the median 
average option. A similar approach to the previous section was 
used for the geometric mean, where however the 2D image was 
build from the log 5/ instead of 5/ to compute the mean value 
(log 5)/. The geometric mean spectrum S gm and the geometric 
standard deviation <x g are given by 



*^em — 



n* 



I IN 



= 10 <log,S/> , cr = 10 V^g^/ASgi 



0] 2 > 



(8) 



The resulting spectra are shown on a log-log scale in Fig. [8] 
for the arithmetic median and in Fig. [9] for the geometric mean 
(extinction-corrected, see below), respectively, except for type 
M where a composite makes no sense as this subsample is too 
small and very heterogeneous. The spectra of different types 
have different normalisations owing to their different wavelength 
coverages. The l<x variation is shown by the shaded area. The 
variation in the flux density across the spectrum is expected to re- 
flect the spectrum-to-spectrum differences caused by differences 
in the continuum shapes and in the properties of the emission 
and/or absorption lines (see Vanden Berk et al. 2001 ). 

The average rest-frame spectra differ significantly from type 
to type, as well as from the Vanden Berk composite of "nor- 
mal" SDSS quasars. The defining spectral features are clearly 
indicated in the median spectra. The median spectrum of the un- 
usual BAL quasars of type A is remarkably dominated by the 
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Fig. 9. Geometric mean composite spectra corrected for intrinsic reddening. The SDSS spectra were combined with JHK fluxes 
form 2MASS (squares with error bars). For comparison, the power-law continuum fit to the geometric mean SDSS composite from 
VandenBerk et al. (2001 ) is shown. 



typic al absorption troughs of FeLoBAL quasars (see Hall et al. 
2002; their table 1), of which the strongest are indicated by the 
dashed vertical lines: Lya, Nv ,11334, Cn ,11334, Siiv ,11398, 
Civ ,11550, Aim ,11860, Fen ,12400, Fen ,12600, Fen ,12750, 
and Mg n ,1280Cq] On the other hand, the spectrum-to-spectrum 
variations for type A are huge. This reflects the particularly high 
degree of diversity among these spectra. As expected, the dom- 
inating spectral features in the median spectrum for type B are 
the absorption troughs from the high ionisation lines O iv ,11033, 
Ly a, N iv ,11240, Si iv ,11398, and C iv ,11550. There are no strin- 
gent indications for BALs in the median spectra of the other 
types, though these features are p resent in the spectra of some 
members of these types (Sect.|4.3|). 



5.6. Intrinsic reddening and corrected absolute magnitudes 

The general trend of BAL quasars to show stronger UV redden- 
ing than non-B AL q uasars and of LoBAL quasars to be even 
redder (see Sect. |5.4| ) is clearly indicated by the type composites 
in Fig. [8] The type F composite indicates significant reddening 
shortwards of Mg II. For type E, we found that quasars with z & 2 
tend to have steeper continua than normal quasars (e.g., Fig. [5]), 
but that the spectra are slightly shallower at higher z. This yields 
an average spectrum for type E with a continuum that reasonably 
fits the SDSS composite of normal quasars. The quasars of the 
types C and D are by definition redder than normal quasars. 

We assumed that the red continua are due to reddening by 
dust related to the quasar and/or the host galaxy and estimated 



the mean intrinsic reddening (E^y) of each type from its geo- 
metric mean spectrum 



12 Most of these lines are actually blends of two or several lines. 
However, in many cases throughout this paper it is adequate to treat 
these blends as single lines. 



(int) > 



logS gm , = logSgm + 0A(E^)[R v + Q x \ 



(9) 



where S gm is the observed geometric mean from Eq. (ISl and 



J gm 

Sgm,o is the emitted geometric mean spectrum before dust ex- 
tinction. As has been shown in several previous studies (e.g., 
Reichard et al. 2003 ), SMC dust provides acceptable fits for the 
intrinsic reddening in quasars. Hence, we adopted the SMC ex- 
tinction curve with R v - 2.93 and Q(A) = E^-v/Eb-v from 
Pei dl992b . We assumed that S gm ,o can be identified with the 
geometric mean spectrum of usual quasars. Vanden Berk et al. 
(2001 ) found that a single power-law adequately fits the contin- 
uum between Ly a and HJ3 with a v = -0.44. We used (E^y) in 
Eq. ^ as a free parameter to fit the continuum of the corrected 
spectrum of each type to the power-law with a v = -0.44 and 
identified the best-fit value with the mean reddening parameter 
of this type. 

This approach is hampered by the UV (A ^ 3000A) con- 
tinua being strongly affected by the unusual absorption and/or 
emission features, perhaps with the exception of type E. On the 
other side, the composite spectra do not cover a wide wave- 
length interval long wards of 3 000 A, except for type F. To stretch 
the useful wavelength interval for the continuum fit, we com- 
bined the SDSS spectra with the fluxes derived from JHK mag- 
nitudes given in the 2MASS catalogue (Skrutzki et al. 2006). 
The percentage of quasars with 2MASS magnitudes varies be- 
tween 45% for type B and 70% for type A with a mean value of 
58%. For each object identified in 2MASS, the JHK fluxes were 
normalised with the same normalisation factors used for the cor- 
responding spectrum. The fluxes derived from the 2MASS mag- 
nitudes generally agree with the extrapolation of the composite 
spectra towards longer wavelengths. 
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The extinction-corrected geometric mean spectra are shown 
in Fig. [9] The mean reddening parameters are given in Tab. [4] 
The values are larger than those found in previous studies 
(e.g., E B - V = 0.08 ( 0.02) for LoBALs (HiBALs) according to 
Reichard et al. 2003 ), which is most likely due to differences be- 
tween the quasar samples. The contributions from host galaxies 
in the optical are obviously suppressed by the strong fluxes from 
the bright quasars. 

As can be seen in Fig.[9| the rest- frame optical colours of the 
unusual quasar types are essentially the same as for the SDSS 
quasar composite (for a similar result see , e.g., Zhang et al. 
2010), while they tend to be redder in the UV. In particular, the 
UV flux of the unusual BAL quasars of type A is substantially 
lower than the normal BAL quasars of type B. This is likely 
caused by quasars for which the UV flux is strongl y sup pressed 
by BAL troughs and/or overlapping troughs (Sect. |6.1| ). Owing 
to these features, the fit to the type A composite is particularly 
uncertain. Types B and F are well-fitted by the SDSS continuum. 
No reasonable fit can be achieved for type D, where the slope 
of the dereddened continuum becomes shallower shortwards of 
Mg 11. Since a substantial fraction of red quasars also show BALs 
(Sect.|4.3|), the combined effect of broad-line absorption may 



have some influence on the average continuum, but there are no 
strong indications of BALs in the type D composite spectrum. 
An extinction curve steeper than the SMC curve would be an- 
other explanation. Finally, the type E geometric mean spectrum 
fits the SDSS composite at longer wavelengths without dered- 
dening but flattens approximately below the C iv line. This may 
be an indication for differences between the weak-line quasars 
of high and low redshift. 

In the next step, we used the mean reddening parameters 
from Tab.|4]to estimate the absolute magnitudes statistically cor- 
rected for intrinsic dust extinction adopting the approximation 



(M (cor) ) 



<Mi>- 



' \ h B-V' 



[Rv 



<2(0>]' 



(10) 



where the angular brackets symbolise the ensemble average and 
ieff ~ ^i,eff/(l + z) i s m e rest- frame wavelength corresponding 
to the effective wavelength of the i band with A^s - 7461 A 
(Schneider et al . 120101) . We again adopted the SMC extinction 
curve from Pei (1992]). As can be seen from the results listed in 
Tab.|4] this correction is significant. The one-sided KS (Sect. |5.2| ) 
applied to the corrected MJ cor) inferred that all six quasar types A 
to F are on average more luminous then the comparison samples 
of normal quasars. The finding of the weaker Civ emission of 
BAL quasars compared to that of non-BAL quasars (Gibson et 
al. 120091) thus appears as a manifestation of the Baldwin-effect 
(Baldwin 119771 Osmer et al. 119941 relation between the C iv 
emission line strength and the UV luminosity. 



5. 7. Intervening absorbers 

Finally, we checked by eye the spectra of all quasars for nar- 
row absorption lines distinct from the quasar lines and at lower 
z than the quasar. One of the easiest species to detect is the 
Mg 11 /L12796,2804 doublet, which is frequently in combination 
with Fe 11 /U2587,2600, Fe 11 ^2374,2382, Fe 11 ,12344, and Mg 1 
/12853. While the origin of the absorbing gas remains unclear, 
a direct association with galaxies is suggested by several stud- 
ies (e.g., Bahcall & Spitzer [196 91 Steid el et al. [19971 Bowen & 
Chelouche 201 1 1 Kacprzak et al. 120 lib . Dust embedded in fore- 
ground absorber systems along the line of sight is expected to 
absorb and scatter the UV photons from the quasar and redden 
the quasar spectrum. Detection of the quasar reddening related to 



absorber systems has been reported e.g. by York et al. (2006) and 
Menard et al. ([2008 ) for Mg 11, by Wild et al. (120061) for Ca 11, and 
by Vladilo et al. (2008 ) for damped Lya absorber systems. While 
the extinction curves of the absorbers are generally considered to 
be consistent with that of the SMC, there are probably also cases 
of Milky Way-like dust extinction that can strongly affect the 
quasar spectrum in a narrower wavelength interval owing to the 
broad extinction bump around 2 175 A (Jiang et al. 2011 ). 

Table 6. Unusual quasars with structures close to the sightline 
on the SDSS images. 



SDSS J 



Object 



AL* 7 Ref. 



024230.65 - 000029.7 


2.486 


B 


NGC 1068 


no 


075518.31 + 245432.3 


1.335 


C 


gal 


no 


081959.80 + 535624.2 


2.235 


C 


gal 


yes 1 


090334.94 + 502819.3 


3.579 


C 


gal group 


yes 2 


103410.26 + 344200.7 


2.503 


B 


? 


yes 


103443.66 + 411120.4 


2.064 


F 


gal 


no 


104513.86+171952.2 


3.679 


E 


gal 


? 


112851.83 + 062315.3 


1.503 


A 


gal 


? 


135048.67 + 433406.2 


1.765 


E 


gal 


no 


143445.95 + 422516.7 


1.407 


F 


gal 


no 


144002.24 + 371058.5 


1.409 


A 


gal 


no 


151150.23 + 140506.4 


2.156 


E 


gal group 


no 


164847.85 + 374006.3 


1.328 


C 


gal 


no 



Notes. (a) Foreground object: gal=galaxy, ?=unclear; ^ Absorption line system 

in SDSS spectrum 

References. (1) Inada et al. f20l0l ; (2) Johnston et al. (2003) 
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Fig. 10. Colour indices g-r (foreground extinction-corrected) as 
function of redshift for all types (dots). Other symbols: quasars 
with indications for foreground systems. 



We identified 76 absorber systems in 66 quasar spectra with 
a mean absorber redshift (z a bs) = 1.44 + 0.45. Among them are 8 
quasars with two and one quasar (SDSS J082747. 14+425241.1) 
with three different absorber systems. The redshift differences 
between the absorber and the quasar, expressed by the velocity 
parameter £ = v/c = (R 2 - l)/(R 2 + 1), where R = (1 + z q )/(l + 
£abs), covers the range J3 = 0.02 to 0.72 (mean value 0.28). The 
fraction of quasars with registered absorption systems strongly 
depends on the spectral type. It is lowest for type A (3%) and 
highest for type E (13%). This is almost certainly a selection 
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bias related to the intrinsic complexity of the type A spectra. 
For any given spectral type, on the other hand, quasars with or 
without these narrow-line absorber systems have on average the 
same peculiarity parameters^ 2 . 



Figure[10| shows the colour index g-r (corrected for Galactic 
extinction) versus z for all unusual quasars (as in Fig. [7}, where 
those showing indications of foreground galaxies are marked by 
symbols. Since absorbers with small J3 may be associated with 
the quasar, we split the sample into two groups of 29 quasars 
with /3 max < 0.2 and 37 with fi max > 0.2. As outflows from 
quasars usually have/3 < 0.2 (Foltz et al. |1983l L the latter group 
(with (J3) = 0.41 ± 0.14) is expected to represent a clean sample 
of quasars with spectroscopically identified foreground galax- 
ies. The majority of our quasars with foreground galaxies are 
redder than the normal quasars from the QCDR7. Hence, we 
cannot exclude that a fraction of the red quasars (types 3 and 4) 
were selected as unusual owing to properties related to the fore- 
ground rather than the quasar itself (as is indeed the case for the 
two gravitationally lensed quasars SDSS J08 1959.80+535624.2 
and SDSS J090334.94+502819.3). However, we did not find any 
significant differences between the reddening properties of these 
quasars in our sample showing evidence of foreground galaxies 
close to the line of sight and those which do not. 

In the Appendix A, we present the results of a simple search 
for extended structures close to the line of sight on direct images. 



6. Very peculiar quasar spectra 

We did not identify any fundamentally new type of unusual 
quasar spectra. However, our selection substantially enlarges the 
sample of known unusual quasars and thus enables statistical 
studies of the various peculiarity types and searches for rela- 
tions between them. Representative examples of our 7 types 
were shown in Fig. [5] The detailed investigation of individual 
spectra is beyond the scope of this paper. However, it is worth 
commenting on a few very peculiar and interesting spectra, in- 
cluding another three mysterious objects and another dozen pos- 
sibly related quasars. 



6.1. Unusual BALs and strong Fe emission 

A selection of SDSS quasars with particularly unusual spectra 
caused by extreme BAL features or/and strong iron emission is 
presented in Fig. jj_ Each panel shows the normalised flux F A (A) 
over the observed wavelength interval 3 900 A to 9000 A, the ob- 
ject name is given at the top of the panel. Here we briefly discuss 
these 12 spectra from top to bottom and left to right. 

The first three panels show quasars with strong 
iron emission. In the first panel, we present SDSS 
J100237.22+270056.£3 which exhibits emission from the 
multiplets Fen ,11785 (UV67,UV191), Fern ,11926 (UV34), 
Fern A2010 (UV48), Fen ,12400 (UV2), Fen ,12600 (UV1), and 
Fen ,12750 (UV62,UV63). The positions of these lines, as well 
as the positions of the typical quasar emission lines (Si iv ,11400, 
Civ ^1542.2,1550.8, Aim ,U1 854.7, 1862.8, Cm] A 1908.7, 
Mgn ^2796.3,2803.5), are marked by the vertical solid lines. 
There is a high-velocity outflow at z a bs = 1.724 seen in Civ, 
Al ii, Al in, and Mg n. The absorption line redward of C iv, close 
to Hen ,11640, is probably a narrow intervening system. SDSS 
J100237. 2+270057 is one of the most luminous quasars in our 
sample with M { = -28.94. 



The spectrum in panel (b) shows a quasar (SDSS 
J09221 1.56+365120.2) where the emission from Fen UV1 and 
UV2 and Fein UV48 is stronger than the (weak) broad emis- 
sion lines Civ and Mgn. Strong iron emission is known to 
be correlated with the occurrence of BAL troughs (Boroson & 
Meyers 1992; Zhang et al. 12010b, as can be seen in many of 
our spectra of types A or F (Sect. 43 ). A nice example is SDSS 
J094225.42+565613.0 in panel (c). The positions of absorp- 
tion troughs at z a bs = 0.735 are marked by the dashed lines. If 
we adopt the systemic redshift of SDSS J094225. 42+5656 13.0 
from the (weak) [O n] ,13728 emission line, the Mg n emission 
line is either blueshifted or absorbed at the red side. In Tab.[3) 
the remark "1st?" is made, indicating possible but not certain 
longwards-of-systemic absorption (see Hall et al. 2002]). The 
Mgn emission of SDSS J094225.42+565613.0 may instead be 
affected in this case by broad Fe n absorption on its red side. 

As a consequence of the overlapping Fe n troughs, the con- 
tinuum of SDSS J094225.42+565613.0 is depressed shortwards 
of Mg ii down to the lowest observed wavelengths. Interestingly, 
the Fe n UV2 trough appears to be deeper than UV1 . As outlined 
by Hall et al. ( 2002 ), the only way of explaining an increase in 
the absorption strength when troughs overlap is spatially distinct 
velocity-dependent partial covering (sdvdpc) of the continuum 
source, provided that the troughs are saturated (as is usually the 
case even when the absorption is not black). Hall et al. empha- 
sise that these objects "may not be all that rare, just difficult to 
recognise", though the only quasar known to them to definitely 
exhibit sdvdpc was FBQS 1408+3054. From the inspection of 
our sample, we found at least another eight candidates probably 
showing that effect (Table [7]). 



Table 7. Candidates for spatially distinct velocity-dependent 
partial covering (sdvdpc) of the continuum source by Fe n UV 
troughs. 



SDSS J 


T 


z 


Remark 


015151.58-093215.3 


A 


1.413 


ot(sdvdpc) 


082350.53 + 244653.1 


C 


1.766 


ot( sdvdpc)? 


094225.42 + 565613.0 


A 


0.831 


ot(sdvdpc?,lst?); sFe(opt) 


101912.84 + 410807.4 


A 


2.450 


nt(sdvdpc);FeinBAL? 


135246.37 + 423923.5 


A 


2.000 


ot(sdvdpc) 


142010.28 + 604722.3 


A 


1.349 


ot( sdvdpc) 


150848.80 + 605551.9 


A 


1.525 


ot(sdvdpc?) 


162527.73 + 093332.8 


A 


1.413 


ot( sdvdpc) 



13 A similar object (not shown here) is SDSS J124244.36+624659.1. 



Notes. Remarks see Table 3. 



In the second row of Fig.[TT] we give three examples of 
LoBAL quasars with possible strong Fein absorption. SDSS 
J021 102.33 - 081007.4 resembles the two possible Fern BAL 
quasars SDSS J014905.28 - 011404.9 and SDSS J081024.75 + 
480615.4 from Hall et al. (2002). It appears to have Fern 
UV34, UV48 troughs (which probably overlap), in addition to 
troughs from Siiv, Civ, Mgn, and weak Fen UV1, UV2 ab- 
sorption. Even more pronounced are the iron absorption troughs 
(plus troughs from other elements) seen in SDSS J101912.84 + 
410807.4. For both spectra, and also for all following quasars 
shown in Fig.[TT] it seems impossible to clearly identify any 
emission line, thus only absorption redshifts can be given. The 
BAL quasars can show, of course, two or more absorption sys- 
tems at different z a bs- The values of z a b s given in Fig. [TT] usually 
refer to the centres of the strongest system. 



18 



H. Meusinger et al.: Kohonen- selected unusual SDSS quasars 



(a) 5D55 J100237. 22+270056. 5 zCem) = 1 . 984 



QJ QJ — 
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Fig. 11. Examples of particularly unusual SDSS quasars: strong UV iron emission (top), Fein absorption (second row), BAL quasars 
with many troughs (third row), and with overlapping troughs (bottom). 



The third row shows another three examples for FeLoBAL 
quasars with many narrow troughs, similar to the three high-z 
quasars of this category discussed by Hall et al. (2002). Although 
the spectra appear to be very complex, a large number of absorp- 
tion features can be identified. Only a selection of the strongest 
troughs is labelled here. We again did not try to estimate the 
emission redshift since emission lines can be mimicked by over- 
lapping absorption troughs that make detailed line identification 
quite difficult. As for the previously known objects of this type, 
no absorption trough reaches zero flux, indicating that there is a 
partial coverage of the continuum source by the absorbers and/or 
a contribution from scattered light (see Hall et al. 120021) . 

The SDSS has discovered several FeLoBAL quasars with 
sharply declining flux near Mgn, which is interpreted as be- 



ing caused by overlapping iron troughs (Hall et al. 12002b . The 
bottom row of Fig. 1 1 exhibits three examples where the ab- 
sorption strength appears to increase from left to right. Here, 
the redshifts refer to the red edges of the absorption troughs. 
Only the longest wavelength lines are marked for the Fe multi- 
plets. SDSS J162527.73+093332.8 displays both many narrow 
troughs at A ^ 7 200 A (observed) and overlapping Fen UV1 and 
UV2 at longer wavelengths (Fen UV2 is deeper than UV1, see 
above). 

SDSS J173049. 10+585059.5 was already described by Hall 
et al. (2002). With the exception of two narrow wavelength in- 
tervals around Cm] A 1908. 7, there is almost no flux shortwards 
of the onset of Fe n UV2 at 7000 A (observed) down to the short- 
est wavelengths. A similar spectrum, but even more extreme, is 
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that of SDSS J0943 17.59+541705.1 at the right-hand side. The 
flux drops rapidly at ~ 9000A (observed) towards shorter wave- 
lengths and is below the detection threshold in the SDSS g and 
u bands, i.e. shortwards of the C iv line. The object is detected 
in 2MASS with K = 14.25. The 2MASS colours correspond 
to a slowly decreasing flux with increasing wavelengths from 
1.2 to 2.2 yum. At first glance, the SDSS spectrum suggests a 
high redshift quasar at z > 6. However, a higher-quality Keck 
spectrum presented by Urrutia et al. (2009 ) clearly indicates the 
Mg ii emission line, along with the typical absorption features of 
a FeLoBAL at z = 2.224. 

Given that SDSS J162527.73+093332.8 is more lumi- 
nous than the other two overlapping trough quasars SDSS 
J173049.10+585059.5 and SDSS J094317.59+541705.1 (M { = 
-26.27 compared to -24.65 and -24.75), all three spectra are 
most likely very similar. The main difference is a slightly higher 
redshift of SDSS J162527.73+093332.8, hence the 3000A jump 
is shifted out of the spectral window. The very red optical-to- 
near infrared colours of SDSS J162527.73+093332.8 (which is 
one of the reddest SDSS quasars; r - K = 1 Al) indicate that 
its flux must strongly increase long wards s of the red edge of the 
SDSS spectrum. Lower-luminosity versions of these quasars at 
z ^ 2.3 are strongly biased against being included in optical 
magnitude-limited samples. 

6.2. Mysterious and possibly related objects 



Hall et al. (120021) discussed at length two objects discov- 
ered by the SDSS with spectra classified by these authors 
as "mysterious": SDSS J010540.75 - 003313.9 and SDSS 
J220445.27+003141.8 (Fig.[l2]b,c). The most prominent spec- 
tral features are: (a) a lack of substantial typical quasar emission 
lines (except broad Fe n emission and [O n] A3130 in the case of 
SDSS J010540.75 - 003313.9); (b) a blue continuum longwards 
of ~ 3200A; (c) a dip at ~ 3000A; (d) a continuum drop-off 
shortwards of Mg n /12800 which appears too steep to be caused 
by dust reddening; and (e) no obvious BAL troughs. Both objects 
show associated Mg n absorption lines, are unresolved FIRST ra- 
dio sources, and are much more luminous than any galaxy. Hall 
et al. discussed various viable explanations for these spectra and 
concluded that none of them are particularly satisfactory. These 
strange objects may be very unusual BAL quasars with partial 
covering of different regions of the continuum source as a func- 
tion of velocity, probably in combination with moderate redden- 
ing. 

Regardless of the physical reason for these mysterious spec- 
tra, they almost certainly represent quasar types that are ex- 
tremely rare in presently available samples. Hall et al. (2002) 
were aware of only two objects with probably similar spec- 
tra, the low-z quasars FBQS 1503+2330 (z = 0.40) and 
FBQS 1055+3124 (z = 0.4 9) from the FIRST Bright Quasar 
Survey (White et al. 120001) . In a search for BLLac objects, 
Plotkin et al. (2008) revealed another two FIRST sources 
with SDSS spectra similar to the two Hall et al. objects 
SDSS J130941.35+1 12540.1 and SDSS J145045.56+461504.2. 
Neither of these two new objects ultimately survived as BLLac 
candidates during the manual inspection by these authors. The 
first similar object that was not targeted as a radio source was 
VPMS J134246.24+284027.5, which was discovered during a 
variability and (zero-) proper motion selection of quasar candi- 
dates (Meusinger et al. 2005 ). All seven of these objects have 
counterparts in the FIRST images with radio fluxes on the mJy 
level and moderate or small radio-loudness parameters. 



Browsing through our sample (Tab.|5]) yields another three 
comparable SDSS objects and about one dozen possibly related 
objects with spectral properties that are partly similar to those 
of the mysterious objects. In addition, we checked a subsample 
of the unknowns (spec_cln=0) in the DR7. Only ~ 10 3 spectra 
of relatively high S/N were clustered in one single SOM. The 
resulting icon map was used solely for selecting spectra that re- 
semble those of the mysterious objects mentioned above. Two 
additional objects of this type were found: the mysterious ob- 
ject SDSS J010540.75-003313.9 from Hall et al. ([2002b and 
SDSS J073816.91+314437.0, which was discussed by Hall et 
al. ([2004)0 

Table[5] lists 9 mysterious objects and another 1 1 objects that 
may represent links between the mysterious and other uncon- 
ventional quasar types ^j The following five objects from this 
list are not included in Tab.[2 VPMS J134246.24+284027.5 
was not targeted by SDSS, SDSS J010540.75 - 003313.9 and 
SDSS J073816.91+314437.0 were not classified as quasars in 
the SDSS DR7, and for both SDSS J130941.35+112540.1 and 
SDSS J161836.09+153313.5 the (wrong) redshifts given in the 
DR7 (z = 4.39 and 4.376, respectively) exceed our upper selec- 
tion limit of z < 4.3 (Sect.07l]). Tab.[S]does not include FBQS 
1503+2330 because the SDSS spectrum does not cover much of 
the wavelength range shortwards of Mgn. The objects that are 
not part of our unusual quasar sample are flagged in the first col- 
umn of Ta b. [8] For clarity, the second column gives the panel in 
Figs. 12 to[l4[in which we plot the spectrum. 

For 1 1 quasars, the redshift could be obtained from emission 
lines. In 9 cases (marked in Tab.[8]), the systemic redshifts were 
computed from the [O n] /13728 line, for the other two spectra 
from Mgn and Cm]. The redshifts of the remaining 9 objects 
were identified with those of the highest-z absorption features. 
For a few objects, e.g., VPMS J134246.24+284027.5, the red- 
shift is fairly uncertain. 

All objects with dubious spec tra w ere checked for proper 
motions (pm) as described in Sect. 4.2 The pm data for VPMS 
J134246.24+284027.5 were taken from Scholz et al. ff997l 
A significant non-zero pm was found for the object SDSS 
J020731. 81+004941.6 (Fig.[3|, which was rejected thereafter 
from the quasar list. The determined pm components [i a cos 6 
and fis are listed in Tab.[5j along with the formal errors and 
the number of epochs N e . The column "pm?" contains our 
conclusion, which is not based on the formal error alone (see 



Sect. |4.2| ). For 13 objects (76%), the results are clearly consis- 
tent with the conclusion that no significant pm is detected. For 
three objects, non-zero pm in one direction cannot be definitely 
excluded. A formally significant pm was measured for SDSS 
J075437. 85+4221 15.3. This result is doubtful, however, because 
the measurements are likely affected by a faint source close to 
the sightline. For the three objects without pm data in Tab.[Sj the 
extragalactic origin is clearly indicated by redshifted emission 
lines. 

The partial spectra in the rest-frame wavelength interval 

The spectra were ar- 
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1350 A - 3800 A are shown in Fig._ 

bitrarily normalised and slightly smoothed by a three-pixel 

boxcar, with the exception of the noisy spectrum of SDSS 



14 Another highly peculiar object from this SOM is SDSS 
J164941. 87+401455.9. This spectrum has a similar type to SDSS 
J033810.85+005617.6 from Hall et al. ([20021 . 

15 Since the selection criteria are not sharply defined, a few more 
quasars may be found with more or less similar peculiar spectra, 
such as SDSS J110711.40+082331.2, SDSS J140800.43+345 124.7, or 
SDSS J100353.68+515457.0. 
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Table 8. Mysterious objects and possibly related quasars. 



Name 



Panel in 
Fig. 12 



Ref. 



Mi 
(mag) 



#i 



E'(intr) 

(mag) 



Ha COS 5 

(mas/yr) 



(mas/yr) 



AL 



pm 



(a) Mysterious objects: 

SDSS J010540.75 - 003313.9 a (c) 

SDSS J085502.20 + 280219.6 (f) 

FBQS J105528.80 + 312411.3 (a) 

SDSS J130941.35 + 112540.1 a (h) 

VPMS Jl 34246.24 + 284027.5 a (i) 

SDSS J145045.56 + 461504.2 (g) 

SDSS J160827.08 + 075811.5 (e) 

SDSS J161836.09 + 153313.5 a (d) 

SDSS J220445.27 + 003141.8 (b) 



1.179 b 

1.511 

0.497 

1.362* 

1.300 

1.894 

1.182* 

1.358* 

1.353 



D 
F 



D 
D 

D 



-26.49 

-25.55 
-25.13 
-26.62 
-25.36 
-26.48 
-27.57 
-26.02 
-27.81 



1.19 

1.16 
0.72 
0.89 
1.17 
1.18 
1.22 
0.49 



0.18 
0.15 
0.45 
0.20 
0.12 
0.30 
0.18 
0.10 
0.22 



8.1+2.1 
-8.1 + 1.0 

-3.1+2.7 
-2.2 + 5.2 
3.7 + 5.8 
3.5 + 5.7 
1.5 + 1.2 
1.4 + 2.1 
-1.6 + 1.3 



-4.1+2.3 

1.5 + 2.0 
-2.2 + 3.5 
-1.2 + 2.6 
-0.5 + 5.9 

1.8 + 1.8 

1.6 + 1.9 
-9.0 + 2.1 
-1.5 + 1.6 



40 
13 

7 
10 

22 



11 

8 

23 






pmx 

pmx 

nopm 

nopm 

nopm 

nopm 

nopm 



id 



8 pmy 
nopm 



(b) Possibly related objects: 

SDSS J073816.91 + 314437.0 a (q) 

SDSS J075437.85 + 4221 15.3 (s) 
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Notes. (a) not in Tab.[3] ^ redshift based on [On]; ^ not in FIRST survey area; ^ "pmx?" and "pmy?" indicate small (~10 mas/yr) proper motion components 
/u a cosS or jug respectively. These results are formally significant but doubtful, as caused by probably erroneous POSS1 positions. (e) "pm?" indicates formally 
significant proper motions in both components that are likely caused by a faint object close to the line of sight. 
References. (1) Hall et al. ( f2002l : (2) Hall et al. ( f2004l : (3) Plotkin et al. (2008 1; (4) Meusinger et al. (2005 1 



J085502.20+280219.6 where a seven-pixel boxcar was ap- 
plied. For comparison, the SDSS quasar composite spectrum 
is shown in the last panel (bottom right). The vertical bars 
above the spectra indicate the typical strong quasar emission 
lines Si iv ^1394,1403, Civ ^1548,1551, Cm] ^1909, and 
Mgn /L12796,2804. The dashed lines below the spectra mark 
the same absorption lines as in Fig.[TT] If the redshift was ob- 
tained from absorption features, the lines above and below the 
spectrum refer to the same z = z a bs- Exceptions are SDSS 
J101723.04+230322.1 and SDSS J073816.91+314437.0, where 
the solid lines mark the highest-z and the dashed ones the 
strongest absorption troughs. 



The three new mysterious objects (Fig. 12 d-f) are SDSS 
J160827.08+075811.5 (z = 1.18), SDSS J161836.09+153313.5 
(z = 1.358), and SDSS J085502.20+280219.6 (z = 1.511). Both 
SDSS J160827.08+075811.5 and SDSS J161836.09+153313.5 
are similar to the Hall et al. object SDSS J010540.75 -0033 13.9. 
A remarkable difference between these two spectra and the 
SDSS J010540.75 - 003313.9 spectrum is the more abrupt tran- 
sition between the blue part and the red part of the continuum. At 
first glance, it is tempting to ascribe the pronounced peak to the 
Mg ii line. However, if we adopt the redshift from the [O n] line 
as the systemic redshift, the peak is clearly placed shortwards 
of 2800A in either case. For SDSS J160827.08+075811.5, the 
redshift of the [O n] line agrees very well with that of an (as- 
sociated?) absorption line system seen in Mg n, Mg i, and Fe n. 
Another new mysterious object is SDSS J085502.20+280219.6, 
which also resembles SDSS J010540.75 - 003313.9 but has a 
stronger depression at shorter wavelengths. 



Similar to the mysterious objects in Fig.[12|a-i, the spectra of 
SDSS J134951.93+382334.1 and SDSS J215950.30+124718.4 
in panels (j) and (k), respectively, show a break in their con- 
tinuum around the position of the Mgn line. The redshift of 



z = 1.516 for SDSS J215950.30+124718.4 was derived from 
a clearly identified narrow Mg n absorption doublet at the posi- 
tion of the break. Both SDSS J134951.93+382334.1 and SDSS 
J2 15950. 30+ 1247 18. 4 are overlapping trough LoBAL quasars 
(see bottom row of Fig.[TT]) and possibly mark a link between 
overlapping-trough objects and mysterious objects. (The dropoff 
at Mgn in SDSS J215950.30+124718.4 is more abrupt than 
in any of the latter objects.) From this point of view, the 
spectra in panels (h) to (k) seem to represent a sequence of 
overlapping-trough absorption from a low covering factor in 
VPMS J134246.24+284027.5 to the highest covering factor in 
SDSS J215950.30+124718.4. 

The spectra in Fig.[T2]m-o show narrower BAL troughs 
with low covering factors. SDSS J151627.40+305219.7 (panel 
p) is very similar to SD SS J073816.91+314437.0 (panel q) 
from Hall et al. (120041 The high-velocity BAL in SDSS 
J101723.04+230322.1 (panel o) makes it plausible that high- 
velocity BALs contribute to the former two spectra. 

Since intrinsic reddening seems to be present in all spectra 
in Fig. 12 we attempted to correct for it and compare the dered- 
dened spectra with the typical quasar spectrum, following Hall et 
al. (2002). As in Sect.|53] we adopted the SMC extinction curve 



at the quasar redshift with i^-v as a ^ ree P arame ter. For z ^ 1.5, 
the object spectrum was fitted to the SDSS composite spectrum 
at wavelengths A ^ 3000A. For the higher-z quasars, we sim- 
ply tried to fit at the longest wavelengths what almost certainly 
amounts to an underestimated amount of reddening. We note that 
the method of deriving E^y is only approximate and the results 
are expected to be rather uncertain, at least for some objects. 
The mean value (E^J) = 0.27 mag of the individually esti- 
mated reddening parameters of the 8 type D quasars from Tab.[8] 
is somewhat larger than the ensemble average of 0.20 mag de- 
rived from the geometric mean composite spectrum of the type 
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Fig. 12. Partial spectra of 9 mysterious objects and 1 1 possibly related quasars. Radio sources detected in FIRST are labeled. Vertical 
bars: positions of typical strong quasar emission (solid) and absorption (dashed) lines, respectively. 
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Fig. 13. As Fig. 12 but for the ratio of the dereddened spectra to the SDSS quasar composite. 
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Fig. 14. As Fig. 12 but for the ratio of the dereddened spectra to the type F composite. 
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D quasars (Tab.|4]and Sect. |5.5| ). For the five type A quasars in 
Tab.[8| we have (E^) = 0.22 mag compared to 0.18 mag for 
the whole sample. 

The ratios of the thus dereddened spectra to the SDSS com- 
posite are shown in Fig.[l3]in the same style as in Fig. [12] The 
individual intrinsic reddening parameters are listed in Tab. [8] 
A perfect fit would result in a horizontal line at the level of 
unity. All dereddened spectra, perhaps with the exception of 
SDSS J09 1940.97+064459.9 in panel (r), show a depression 
of the flux relative to the composite at short wavelengths. In 
most but not all cases, this sharp decline sets in at wavelengths 
around the Mgn line. This behaviour may be caused by very 
wide iron absorption troughs that overlap at substantially dif- 
ferent velocities (Hall et al.|2002). This depression is also seen 
in the spectra in rows 4 and 5, where the presence of absorp- 
tion troughs from Mgn and Fen UV1 and UV2 (plus Sin and 
Fern UV3£k UV48 in the case of SDSS J120337.91 + 153006.6, 
panel (n)Jjis more obvious. The combination of reddening and 
partial-covering overlapping-troughs could perhaps explain such 
objects as SDSS J145045.56+461504.2 (panel g) and SDSS 
J010540.75 - 003313.9 (panel c), while narrower troughs could 
perhaps explain the spectra in panels (m) to (o). As a remarkable 
exception, the unusual spectrum of SDSS J09 1940.97+064459.9 
can be explained by the 2 175 A absorption trough at the quasar 
redshift (at which narrow associated Mg n absorption is seen). 

Hall et al. (2002) discussed, among others, the possibility 
that the strange continua could represent reddened versions of 
quasars with weak Mgn but strong broad Fen emission. To 
check this idea for the quasars in Fig. 12 we repeated the dered- 
dening procedure with two modifications. First, the composite 
spectrum of strong iron emitters (type F) was used as refer- 
ence instead of the SDSS quasar composite. Second, the indi- 
vidual quasar spectrum was fitted to the reference spectrum at 
either end. Fig. 14] shows the ratios of the dereddened spectra 
to the reference spectrum. A good fit was achieved for SDSS 
J091940.97+064459.9. However, for most of the spectra the de- 
pression in the observed spectrum shortwards of Mg n was ob- 
viously not eliminated. Moreover, strong dereddening produces 
an unexplained bump at 2700A. As argued by Hall et al. (2002), 
Fe ii emission may be present longwards of Mg n (near 3200A) 
but is not expected immediately shortwards of Mg n. In agree- 
ment with Hall et al., we conclude that the majority of the spectra 
cannot be explained by reddened strong iron emission alone. 

the shape of the 



For several objects in Figs. 12 to 14 



spectrum around Mgn resembles a broad double-line profile. 
The most interesting case is SDSS J101723.04+230322.1. Very 
broad double- shouldered Mgn emission, in combination with 
reddening and broad, shallow absorption troughs, is another pos- 
sible scenario mentioned by Hall et al. (2002) to explain their 
two mysterious objects. As can be seen in Figs.[l3]and[T4] the 
spectrum of SDSS J101723.04+230322.1 cannot be understood 
as being caused by these double- shouldered line profiles alone 
or in combination with dust reddening, strong iron emission 
and associated narrow-line absorption. In addition, overlapping 
broad absorption troughs from various elements, in particular 
Fe ii, are definitely needed to explain the complex shape of this 
spectrum. Remarkably, the apparent double-line profile in SDSS 
J101723.04+230322.1 is seen not only in the Mgn line but also 
for Cm] and Civ. Only a few quasars have been reported so 
far to definitely display this structure in other lines than Ha and 



16 SDSS J152438.79+415543.0 and SDSS J120337.91 + 153006.6 are 
anoth er t wo sdvdpc-candidates with deeper Fen UV2 than UV1 
(SectlST 



H0 (Halpern et al. 2001 ; Strateva et al. [20031 Luo et al. [20091 
Chornock et al. 120101) . The selection of double-peaked AGNs 
has been primarily based on B aimer lines (Halpern et al. 120011 
Eracleous et al. 120041) . Ultraviolet observations have shown that 
their high-ionisation lines (e.g., Civ) frequently lack double- 
peaked profiles. However, double-peaked AGNs appear to be 
rather inhomogeneous as a class, thus we feel that this difference 
does not conclusively rule out the double- shouldered hypothesis 
for SDSS J101723.04+230322.1. It would be interesting to study 
the line profiles of the B aimer lines via infrared spectroscopy. 

Most of the quasars in Tab.[8]are FIRST sources, but radio- 
weak, i.e. close to the threshold between radio-quiet and radio- 
loud quasars. The mean radio loudness of the seven radio- 
detected mysterious objects is (R{) = 0.97 + 0.28. Among the 
three new mysterious objects, only SDSS J160827.08+075811.5 
was detected by FIRST. SDSS J085502.20+280219.6 is 2.6 mag 
fainter in i than SDSS J160827.08+075811.5 (and a FIRST flux 
as small as 0.2 mJy thus still corresponds to R x = 0.4). The third 
quasar, SDSS J215950.30+124718.4, is not in the FIRST area. 
The mean radio loudness parameter of the 15 FIRST sources 
in Tab.[8]is (R{) = 1.02 (+0.41). The strongest radio source, 
SDSS J091940.97+064459.9, was classified by Plotkin et al. 
(2008) as a lower-confidence BL Lac object. None of the objects 
from Tab.[8] could be identified in the ROSAT All-Sky Survey 
Faint Source Catalog (Voges et al. 2000). Weak X-ray emission 
is a characteristic property of BAL quasars and is usually at- 
tributed to strong X-ray absorption, which enables and enhances 
the acceleration of the BAL outflow driven by radiation pres- 
sure ( Murray & Chiang 1998; Proga et al. 120001 G allagher et al. 
120061 Gibson et al. [20091 Streblyanska et al. l20T0l) . 

7. Summary and conclusions 

We have compiled a catalogue of unusual quasars from the un- 
precedented spectroscopic database of the SDSS DR7. Unusual 
quasars were selected from ~ 10 5 spectra with z = 0.6 to 4.3 
classified as quasars by the SDSS spectroscopic pipeline. The 
selection method is essentially a combination of the power of 
Kohonen 's (Kohonen 1982 2001 ) self-organising maps and the 
detailed visual inspection of the selected spectra. We paid partic- 
ular attention to reject contaminants, i.e. spectra from rare spec- 
tral types of other astrophysical objects (e.g., white dwarfs and 
hybrid spectra of two objects within the fibre aperture). Proper 
motions were checked first by cross-correlating our sample with 
the PPMXL catalogue (Roser et al. 120101) . which contains ~ 90% 
of the selected objects. In addition, and more importantly, we 
used all available multi-epoch positions from different sky sur- 
veys to estimate the proper motions of 121 selected sources. The 
results of the pm determination led to the rejection of seven ob- 
jects from the quasar sample. 

The final catalogue contains 1005 quasars, which were clas- 
sified into 7 different types: (A) LoBALs and unusual BALs 
(21%), (B) strong normal BALs, HiBALs (21%), (C) red quasars 
(12%), (D) quasars which appear red shortwards of ~ 3000 A 
only (15%), (E) weak-line quasars (18%), (F) strong iron emit- 
ters (11%), and (M) miscellaneous. Our approach was primarily 
aimed at striking outliers and the catalogue is expected to be 
largely complete in this respect. However, the selection is not 
based on sharply defined quantitative criteria and the selected 
sample is thus incomplete in a quantifiable sense. Nevertheless, 
the compilation is expected to be very useful for studying rela- 
tions between the various types of spectral peculiarities or se- 
lecting particularly interesting individual objects for detailed in- 
vestigations. 
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The analysis of this sample yields the following conclusions 
for the main types A to F: 



- The arithmetic median composite spectra (Sect. |5.5| ) clearly 
differ from type to type and from the SDSS composite spec- 
trum of normal quasars. For type A, the huge spectrum-to- 
spectrum variation reflects a high level of diversity. The spec- 
tral peculiarity, measured by the mean square deviation of 
the individual spectrum from the SDSS composite, is on av- 
erage largest for types C,D, and A and smallest for type E. 
Unusual BALs are frequently accompanied by red continua 
and strong iron emission. 

- We combined the SDSS spectra with the fluxes from 2MASS 
to compute geometric mean spectra with significantly im- 
proved wavelength coverage (Sect. |5.5| ). These compos- 
ite s pectr a were used to estimate the intrinsic reddening 
(Sect. |5.6| ). In agreement with previous studies, we found 
that BAL quasars have, on average, significantly redder con- 
tinua than normal quasars with E(B - V) ~ 0.1 mag for 
type B adopting SMC extinction. The unusual BAL quasars 
(type A), which are mostly LoBALs, are even redder with 
E(B - V) « 0.18 mag. The decline in the dereddened type 
A composite in the UV is ascribed to the combined effect of 
the UV absorption troughs. The strongest mean reddening is 
found for the type C with E(B - V) ~ 0.38 mag. On the other 
side, weak-line quasars are on average insignificantly redder 
than normal (E(B - V) « 0). 

- For the quasars with substantial reddening in the UV only 
(type D), the continuum slope of the SDSS quasar composite 
could be fitted only with an extin ction curve that is steeper 
in the UV than for the SMC (Sect.[5^|). 

- After the (statistical) correction for intrinsic extinction, the 
KS test for the absolute magnitudes yielded that all six types 
are on average more l umin ou s tha n comparison samples of 
normal quasars (Sects. |5^2| and |5.6| ). 

- The fraction of quasars with radio detections increases with 
the spectral peculiarity. This reflects most likely a selection 
bias: a substantial fraction of the quasars with very unusual 
spectra were selected for foloow-up spectroscopy because 
they are radio sources. On the other hand, unusual BALs 
(as well as quasars with strong iron emission) tend not to 
have high radio luminosities. A similar, though not statisti- 
cally significant trend was found for the red quasars. It is thus 
tempting to speculate that there exists a larger population of 
unusual quasars that have not yet been disc overe d simply be- 
cause they are too faint radio sources (Sect. [53}. 

- To estimate the impact of foreground galaxies along the line 
of sight on the spectral peculiarities, we have checked both 
the SDSS images and spec tra for finger prints of such in- 
tervening matter (Sect. |5.7| ). About 1% of the quasar im- 
ages show nearby (usually faint) structures and for ~ 8% 
we detected narrow absorption lines that are distinct from 
the quasar lines. There are at least two quasars in the sam- 
ple where the spectral peculiarity is clearly caused by fore- 
ground galaxies. Moreover, two of nine high-z quasars iden- 
tified on images from the Hubble Space Telescope show ex- 
tended structures very close to the sightline (Appendix A). 
Nevertheless, we did not find any conclusive evidence that 
the sample properties are significantly affected by the extra- 
galactic foreground. 

- Our sample contains several objects with very peculiar spec- 
tra, among them FeLoBAL quasars with large numbers of 
narrow absorption troughs or extremely wide overlapping 
troughs. The most spectacular example of the latter type is 



SDSS J0943 17.59+541705.1 at z = 2.22 where the contin- 
uum is strongly depressed from Mg n at red edge of the spec- 
trum all the way down to Lya at the blue edge (Sect. |6.1| ). 
- We constructed a small sample of nine quasars with spectral 
properties similar to the two "mysterious" objects discov- 
ered by Hall et al. ('2002') and another 1 1 more or less similar 
quasars (Sect. 62). The majority (75%) of these objects are 
radio sources, but radio- weak (R{ « 1). Both moderate red- 
dening, FeLoBAL features, and strong Fe emission were fre- 
quently indicated. This combination may explain the shape 
of the continuum with the characteristic drop-off shortwards 
of the position of the Mg n line but does not account for the 
lack of typical quasar emission lines and the shape of the 
spectrum around the position of the Mg n line. Broad double- 
peaked emission, in addition to other effects, is an attractive 
idea (Hall et al. 2002). Double- shouldered Mgn emission 
seems to be indicated in some spectra. The best case is SDSS 
J101723. 04+230322.1 where similar profiles are, however, 
seen as well around Civ and Cm], which is quite unusual 
for the known double-peaked quasars. Moreover, we have 
also presented a few quasars where the continuum drop-off 
appears at much shorter wavelengths, which seems difficult 
to explain in this way. A satisfactory explanation of these 
puzzling spectra remains a challenge. 

Acknowledgements. We thank the referee, Patrick B. Hall, for providing 
many constructive comments and suggestions that significantly improved this 
paper. This research is based on the Sloan Digital Sky Survey (SDSS). 
Funding for the SDSS and SDSS-II has been provided by the Alfred P. 
Sloan Foundation, the Participating Institutions (see below), the National 
Science Foundation, the National Aeronautics and Space Administration, 
the U.S. Department of Energy, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Council for England. 
The SDSS Web site is http://www.sdss.org/ The SDSS is managed by the 
Astrophysical Research Consortium (ARC) for the Participating Institutions. 
The Participating Institutions are: the American Museum of Natural History, 
Astrophysical Institute Potsdam, University of Basel, University of Cambridge 
(Cambridge University), Case Western Reserve University, the University of 
Chicago, the Fermi National Accelerator Laboratory (Fermilab), the Institute for 
Advanced Study, the Japan Participation Group, the Johns Hopkins University, 
the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle 
Astrophysics and Cosmology, the Korean Scientist Group, the Los Alamos 
National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max- 
Planck-Institute for Astrophysics (MPA), the New Mexico State University, the 
Ohio State University, the University of Pittsburgh, University of Portsmouth, 
Princeton University, the United States Naval Observatory, and the University of 
Washington. 

This publication makes use of data products from the Wide-field Infrared 
Survey Explorer, which is a joint project of the University of California, Los 
Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, 
funded by the National Aeronautics and Space Administration. This research is 
based also on observations with the NASA/ESA Hubble Space Telescope, and 
obtained from the Hubble Legacy Archive, which is a collaboration between 
the Space Telescope Science Institute (STScI/NASA), the Space Telescope 
European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy 
Data Centre (CADC/NRC/CSA). Finally, this research has made use of the 
SIMBAD database, operated at CDS, Strasbourg, France. 



Appendix A: Images of unusual quasars in the 
SDSS and in the Hubble Legacy Archive 

Unusual spectral properties can be the result of a positional co- 
incidence with a foreground object in combination with gravi- 
tational lensing (Irwin et al. 119981 Chartas l2000l) . We used the 
SDSS Explorer and found that for ~ 20% of the quasars, with 
only slight fluctuations from type to type, additional objects are 
seen within 12"x 12" fields centred on the quasar. A substan- 
tial fraction is expected to be faint stars from our galaxy. For 
1.3% of the quasars, we found another object closer than ~ 2" 
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or extended structures that likely overlap the line of sight to- 
wards the quasar. In this case, the fraction is largest for type C 
(3.4%), but the samples are much too small for a statistical anal- 
ysis. These 13 quasars are listed in Tab.[6] Among them are the 
two known lensed quasars SDSS J08 1959.80+535624.2 (Inada 
et al. l20T0t and SDSS J090334.94+502819.3 (Johnston et al. 
120031) . The SDSS image of SDSS J081959.80+535624.2 shows 
a double structure with a separation ^ 4". The spectrum indi- 
cates that there is an absorption line system at z a bs = 0.294. 
SDSS J090334.94+502819.3 was classified as an unusual quasar 
of type C because of its red continuum between Ly a and C iv. 
In addition, the C iv emission line is absorbed at the red wing. 
However, as demonstrated by Johnston et al. (2003), decom- 
posing the spectrum into a luminous red galaxy component at 
z = 0.388 and the background quasar component yields a normal 
quasar spectrum. Another interesting case is the quasar SDSS 
J024230.65 - 000029.7, which is seen through an outer spiral 
armofNGC1068. 



Table A.l. Unusual quasars on 


HSJ 


'images. 




SDSS J 


z 


T 


Structure 


Ref. 


005006.34-005319.3 


4.347 


B 






090334.94 + 502819.3 


3.579 


C 


gal. group 


1 


114756.00-025023.4 


2.560 


C 


faint nearby 


2 


122622.03 + 662018.0 


3.874 


B 






134026.43 + 634433.1 


2.784 


A 


faint 5" 




150424.98 + 102939.1 


1.835 


E 


gal ~ 4" 




155633.78 + 351757.3 


1.499 


A 






173049.10 + 585059.5 


2.033 


A 






220557.02 + 121239.7 


1.382 


E 







References. (1) Johnston et al. f2003l ; (2) Bentz et al. J20081 



One of us (B.K.) has performed a systematic search of the 
Hubble Legacy Archive for images from the Hubble Space 
Telescope covering q uasar s from our catalogue. Nine detec- 
tions were found (Tab. |A.l| ). Two quas ars c learly show extended 
structures close to the sightline (Fig. |A.l| ). The gravitationally 
lensed quasar SDSS J090334.9+502819 mentioned above is 
seen through the inner region of a luminous featureless galaxy. 
The other object, SDSS Jl 14756.0 - 025023 at z = 2.560, is 
connected with a faint structure, which may be at the redshift 
of the quasar (Bentz et al. 2008). Both spectra show red con- 
tinua of type C and line absorption at z a bs ~ z em - Extended struc- 
tures are also seen around SDSS J134026.4+634433 and SDSS 
J150425. 0+102939, although at somewhat larger distances of 
~ 4" from the line of sight. There are no clear cases of the sig- 
natures of foreground objects in the spectra of these quasars. 



Note added in proof. After this paper has been accepted for 
publication, we became aware of the paper by Plotkin et 
al. (120101) on optically selected BL Lac candidates from the 
SDSS. As a side-product, these authors discuss a small sam- 
ple of serendipitously recovered/discovered higher-redshift ob- 
jects that show extreme drop-off s in their continua bluewards of 
restframe 2800 A. We realised that eight out of their nine ob- 
jects are part of our sample of mysterious and related objects. 
Among them are four objects (SDSS J160827.08+075811.5, 
SDSS J161836.09+153313.5, SDSS J134951.93+382334.1, and 
SDSS J215950.30+124718.4) which were supposed to be newly 
discovered by the present study. 



•A 



(a) SDSSJ090334.9+502819 z = 3.579 



(c) SDSS JB4026.4+634433 z = 2.784 



'^N& 



(b) SDSS JL 14756.0-025023 z = 2.560 



♦ 



(d) SDSS JL50425.0+102939 z = 1.835 



Fig. A.l. 10" x 10" cutouts from the Hubble Space Telescope 
images of four quasars (unresolved objects close to the centre) 
showing extended structures within 5" from the line of sight, ((a) 
to (c): ACS, (d): WFC3). 
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